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1.0  INTRODUCTION 


V  .  .  , 

■  The  work  presented  in  this  report  is  a  continuation  of  work  started 
under  NASC  Contract  No.  00019-68-C-0274,  Collocation  Flutter  Analysis  Study. 
This  wOrk,  which  is  presented  in  three  volumes,  is  to  update  and  document 
computer  automated  flutter  analysis  techniques.  The-volumes- -contain:  ’  ^  / 


Volume  II 


Volume  III 


The  report  contains  a  set  of  instruction  manuals  with  sufficient  infor¬ 
mation  to  operate  each  program.  The  programs  are  coded  in  Fortran  IV,  and  require 
only  a  minimum  amount  of  modification  to  be  operable  on- most  computers ; 

The  two  programs  presented  in  Volume  I,  the. Subsonic  Strip  Theory 
Unsteady  Aerodynamics  Program  and  Supersonic  Piston  Theory  Unsteady  Aerodynamics 
Program,  calculate  unsteady  aerodynamic  influence  coefficients,  AICs;  These 
AICs  can  be  used  directly  as  input  for  the  Collocation  Flutter  Analysis  Program 
presented  in  Volume  IV  of  Reference  1;  and  when  transformed  into  generalized 
aerodynamic  forces,  they  can  satisfy  the  input  requirements  for  the  Modal  Flutter 
Analysis  Program  described  in  Volume  III  of/  this  report.  The  programs  were 
developed  using  the  strip  theory  approach  so  that  chordwise  camber  could  be  in¬ 
corporated  into  the  analysis.  The  subsonic  program  is  applicable  to  wings  of 
moderate  to  high  aspect  ratios.  The  supersonic  program  is  applicable  to  wifgs 
of  all  aspect  ratios;  however,  for  wings  of  low  aspect  ratio,  it  is  recommended 
that  the  computer  program  presented  in  Reference  2  be  used.  This  program, 
uses  a  normal  mode  analysis  technique;  which  is  inherently  less  accurate 
than  the  AlC-Colldcatiqn  Method.  However,  the  deformation  modes  of  wings  of  low 
aspect  ratio  are  more  accurately  described  in  the  modal  analysis.  It  is  for  this 
reason  that  the  normal  mode  method  is  recommended  for  wings  of  low-aspect  ratio. 


-  ^ubsonic  Strip  Theory  Unsteady  Aerodynamics  Program 

(Strip).' 

; ^Supersonic  Pisfcoh  Theory  Unsteady  Aerodynamics  Program 
(Piston) „  ^ 

-  Unsteady  Aerodynamics  Generalized  Force  Programs  for 
the  Subsonic,  Sonic  and  Supersonic  Flight  Regimes 

-  Structural  Analysis  Program  -  FLUENC-100C 
Component  Mode  Synthesis  Program  (COMSYN) 

Modal  Flutter  Analysis  Program  (MOFA) 


Volume  I 


(i) 


"i—  - 


2.0  SUBSONIC  STRIP  THEORY  AERODYNAMICS  PROGRAM 

2.1  THEORETICAL  DEVELOPMENT 

It  is  desirable  to  consider  the  derivations  of  the  AIC's  in 
the  simplest  form,  that  is  from,  a  strip  theory  approach  in  which 
the  flow  along  any  section  of  the  wing  can  be  considered  two  dimen¬ 
sional;  The  derivation  for.  the  rigid  chord  shown  is  taken  from 
Reference  3.  The  derivation  for  the  flexible  chord  is  an  extension 
of  the  rigid  chord  case  and  is  developed  in  ,a  parallel  manner.  Three 
basic  relationships  must  be  established  to  obtain  AIC’s  from  any  aero¬ 
dynamic  theory;  they  are  (1)  the  pressure-downwash  relation;  (2)  force- 
pressure  -.elation;  arid  (3)  downwash-deflection  relation.  For  the  strip 
theory  ci se  the  pressure-downwash  relationships  are  available  in  an 
equivalent  form.  Theodorsen  (Ref.  4)  has  integrated  the  pressure 
relationships  and  has  presented  the  above  information  as  a  tabulation 
of  oscillatory  coefficients;  L,  M,  N,  T.  The  force-pressure  relation¬ 
ship  is  established  through  the  oscillatory  coefficients  and  a  corelation 
of  the  force  systems  in  Figs.  2.1.1  a  &  b.  The  downwash-deflection 
relationship  is  established  through  the  geometrical  relationships  shown 
In  Figs.  2.1.1  a  &  b. 

In  the  case  of  strip  theory,  the  matrix  of  AIC's  appears  in  a 
partitioned  form;  For  example,  the  AIC's  for  the  two-strip  wing  appear 
as 


St  - 


where  the  C^.  are  the  AIC's  for  strip  "i”.  Thus  it  is  only  necessary 
to  derive  ingeneral  form  the  AIC's  for  one  strip.  This  is  then  applied 
to  each  strip,  and  the  complete  matrix  is  compiled  as  shown  above. 

A  survey  of  two-dimensional  oscillatory  aerodynamic  theory  yielded  the 
incompressible  solutions  of  Theodorsen,  Ref.  4,  and  of  Theodorsen  and 
Garrick,  Refs.  5,6,  and  the  tabulations  of  Smilg  and  Wasserman,  Ref.  7. 
The  subsonic  solution  has  been  tabulated  by  Timman,  Van  de  Vooren,  and 
Griedanus,  Ref.  8.  These  solutions  are  for  the  case  of  a  rigid  airfoil 
and  control  surface.  The  incompressible  solution  for  the  case  of  a 
flexible  chord  undergoing  parabolic  changes  in  camber  has  been  obtained 
by  Spielberg,  Ref.  9.  The  incompressible  case  for  the  cambering  airfoil 
with  control  surface  was  solved  by  Tyler,  Ref.  10. 


Rigid  Chord 

The  theory  is  developed  for  the  general  case  with  control  surface,  and 
degenerates  to  a  second  order  matrix  (upper  left  partition)  for  the  non- 
contri 1  surface  configuration. 


DEFINITION 


The  force-pressure  relationship  derived  from  Fig.  2.1.1  is 

r,  ,  i  i 


1  (d+c  ), 
j  a 


2.1.2 


from  Ref.  4,  the  oscillatory  coefficients  are  defined  as 

iLi  f1  “I'm  «■„  i si 


(M  (  =  ?r  cos  Ap  u^b^'Ay 


1°  V  Ma 


0  b 


2.1.3 


The  geometrical  relationship  between  the  downwash  and  the  deflection  as 
derived  from  Fig.  2.1.1  is 


hi  1  0  •  o  /  h.  \ 

'  1  1  J 

ba  =  -b/d  b/d  j  0  jh^.t- 

bg  b/d  -(b/d+b/c  )  b/c  (h_| 

cL  3.  \  O  I 


2.1.4 


L  Substituting  2.1.4  into  2.1.3  and  the  result  into  2. 1.2, inverting  and  multiplying 
we  get 


F2/  =  17  cos^ 1 


1 _ 

-b/d  | 

i 

0 

b/d  j 

.0 

1 

1 

1  o 

1 

1b  ‘-i1- 


T.  T  tf  J  _b/d  -(b/d+b/c  ) 
nap  3 


B/c  J 
-  a 


2.1.5 


% 


Comparing  2.1.5  with  2.1.1  we  see  that 


ir  cos  A  (b/br)2(Ay/s)  0  b/d  1 


-*/d  |  b'd  1  ^  La  !  Vjr  3 

!  ! 

b/d  I  -(b/d+b/c  )  K  H  ill  -b/d 
j  a  n  a  j  B 


0  b/c 


lTh  Ta  W/d 


'(b/dfb/c  )  b/c 


Flexible  Chord  with  Parabolic  Camber 


2.1.6 


The  theory  Is  developed  for  the  general  case  with  control  surface  and 
degenerates  to  a  third  or4er  matrix  (upper  left  partition)  when  the  control 
surface  is  absent. 

Using  virtual  work  to  establish  the  force-pressure  relation 
SW  =  l$hc/4  +  M6a  +  N6;  +  16  f 


2.1.7a 


6hcM\T 


b6a 


j  b6B 


2.1.7b 


6W  =  FJh.  +  F06h0  +  F'6h_  +F.6h. 

11  22  33  4  4 


2.1.8a 


6hi  T 


i-2 


K 
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2.1.8b 


Substituting  2.1.11b  into  2.1.10  to  obtain 
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therefore 

[%]  -  *  aZaA?*)  [a]  1  [t]  |a 


Development  of  the  [a]  matrix. 


1 


Figure  2.1.2 


Using  Lagrangian  Curve  Fit  with  Figure  2.1;'2 


hx  (x-x2):  (x-x3) 
(x1-x2)  (x1-x3) 


0^3)  +  „3 

(x2-x1)  (x2~x^)  j  (x3~x1)  (x3~x2) 


2.1.15 


n  (x) 


2x  -x2  ~x3 
ll  (x^)  Cx1-x3) 


,  .  2x-x,-x_ 

•  “2^  v  -  . _ 

<*27xz>  (x2"x3) 


2x-x^-x2 
(x3“x1)  (x^-SQ 


2.1.16 


We  can  now  relate  1^,  h2,  h3>  to  h^,  ba,  c  ,  b{3 


h£e~  {h£e/hl}  hl  +  {h£e/h2}  h2  *  {h£e/h 3*  h3 


hc/2  {|lc/2/hl}  hl  +  {hc/2/h2}  h2  +  {^c/2/h3}  h3 


hte-  =  {hte/hl}  hl  +  <  W  h2,  *  {hte/h3}  h3 


2.1.17 


where  we  have  used  the  following  notations 


<  W  -  (V-x2)  Cx1  -x3) 

1  Vh2>  =  Ule-\)  <V*3>  (x2  "V  (x2  "V 

lW  ■  (x2e'xl)  <x3  "V  (x3  ‘V 


{hc/2/hl1 

v-< 

it 

(xc/2-’b) 
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-*3) 

!hc/2^h2^ 

■  <’Ic/2-xl) 

(ltc/2-x3>  (x2 

-3^)  (x2 

“X3) 

*hc/2^h3^ 

"  (xc/2-xL' 

Cxc/2'XL)  (X3 

-xp  (x3 

-X2) 

<hta/hl) 

■  (VV 

^xte~^  CxL  ' 
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V 
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■ 

‘W  tX2  - 
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V 
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- 
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X2) 

2.1.18 


(8) 

f- 
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E.-'-k  - 


hc/4  "  14  (3\e  +  hte> 


“  {hc/4/hl*  \  +  ^c/4^h2}  ^2  +  {hc/4/h3?  h3 


2.1.19 


where 


hc/4/hi  “(3/4)h£e/hi  +  a/4)  hte/hi*  1  “  1>2»3 


ba-  *  d/2)  Chte-h£e) 


{Bo/h^lh,  +{ba/h2  }h2  +  {ba/h^}  hg. 


where  bcs/h^  *  (1/2)  (h^/t^  -h£e/ht)»  i  “  1,2,3 


2.1.20 


2.1.21 


2.1.22 


«  "  hc/2  '  (1/2)  (h  e  +  hte> 


=  fe  /h±}  hx  +{t/h2l  h2  '+  (c/h3}  h3 


2.1.23 


where 


5/^  =  hc/2/hi  ~  ^hte^l  +  hte^hP  *  1  1,2,3 


2.1.24 


.  x„  -x„  X  -  X, 

.  ' ?  ,x  x  3  2 _ .  3  1  . 

“  4  3;  =  (xa  -x2)(Xl  -  x.3)  ni  +  (x2  -  x1)(x2  V  x3)c  n2 


.  2x  -  x-  -  x0  , 
rt*  3  1  /  n. 

Cx3  -Xx)  (x3  -x2)  - 


2.1.25 


bB  =  b(h^  -h3)  +  bh'  (x3) 


2.1.26 


b(-*  3  ~  *2^  .  hx  +  b(~x3  ~xi)  h2  +  b(2x3  ~X1  -  1  h  +  b  h4 

(xl-x2'*  ^xl-x3^  (x2  ~x^).(x2~x3)  (x3  -xx)  Cx3  -x2)  ca  3  c& 


i  »  b{h'c»3)/h1)  h3  +  b  {h  (»3>/h.)  +  b  {h  (x3)/n3)  h3  + 


<»'  f»3)/h,]  -  (*3  -  V  .  ■■  - 

lb  (xj/fij  -  fe3  "M? _ 

(xL  -x^)  (X2  -x^) 


{h  (x3)/h3>  «  2x3, ~  *x  ~  x2 

(x3  -  xi)(x3  -x2) 


{h  (x3)/h3>  =  1 


•  o|o* 


Thus-  we  may  write  the  geometrical  relation  matrix. 


l  ^ 

fhc/A! 

{hc/4^h2} 

{1W*/h2} 

0 

j  ba 
\ 

l  =  ’ 

f 

b{a/h^}’ 

b{a/h2} 

b{a/h3) 

0 

h2 

\ 

)  { 

> 

{?/h2} 

(C/h3} 

0 

h3 

,b6 

r 

b(h'  (x3/h1?  b(h*  Gc3)/h2)  b{h*  (x^/hg} 

b{h’  6«3)/h^} 

n 

where  the  elements  of  [A]  are  defined  in  equation  2=1*29 


(ID 


2.2 


PROGRAM  DESCRIPTION 


A  general  program  to  calculate  a  set  of  aerodynamic  influence  co¬ 
efficients  using  incompressible  strip  theory  has  been  developed.  The  method 
is  applicable  to  wings  of  moderate  to  high  aspect  ratio  and  speeds  in  the 
subsonic  regime.  The  analysis  can  be  performed  for  wings  with  a  rigid  chord 
or  a  flexible  chord.  The  effects  of  a  flexible,  chord  are  accounted  for  by 
the  introduction  of  parabolic  cambering  if  the  bending  inode  is  parabolic 
and  the  torsion  mode  linear  in  the  region  surrounding  the  strip  under  con¬ 
sideration.  The  analysis  can  be  performed  with  or  without  a  control  surface. 
The  method  used  is  based  upon  the  most  fundamental  solution  in  unsteady  flow 
by  Theodorsen  for  the  oscillating  two-dimensional  airfoil  in  an  incompressible 
flow  and  the  extensions  to  include  damber  by  Speilberg  and  Tyler.  The  steady 
state  case  is  available  as  a  limiting  case  of  the  oscillating  case  for  use  In 
static  aeroelastic  analysis.  The  AICs  relate  the  aerodynamic  forces  to  the 
surface  deflections  through  the  following  definitions.  In  the  oscillatory 
case. 


(?) 


p  ?*2br2  Ch]  { h} 


and  in  the  steady  case,. 


(F> 

s 


=  (1/2)  pv2 


The  AICs  are  derived  for  each  strip  considering  the  airfoil  to  have 
up  to  four  degrees  of  freedom:  pitching,  plunging,  cambering,  and  control 
surface  rotation.  The  program  provides  the  AICs  in  printed  and  optional 
punched-card,  output  format.  The  punched-card  output  satisfies  the  input 
requirements  cf  the  Collocation  f  lutter  Analysis  Program  (Ref.  1) .  The 
program  capacity  is  25  surface  strips  and  50  values  of  reduced  velocity; 


2.2.1  PROCESSING  INFORMATION 


A.  OPERATION 

Standard  FORTRAN  IV  processor  system  .  Operable  on>  the  GE  635 
computer. 

B.  CORE  STORAGE 

The  program  STRIP  requires  a  minimum  of  20,000  memory  units  for 
execution. 

C.  ADDITIONAL  MACHINE  COMPONENTS 

Standard  FORTRAN  input  tape  (5) 

Standard  FORTRAN  output  rrint  tape  (6) 

Standard  FORTRAN  output  punch  tape 


(12) 


* 


2.3  INPUT  INSTRUCTIONS 
UNITS 

Since  all  of  the  input  dimensions  are  geometrical  and  the  aerodynamic 
matrix  is  dimensionless,  only  a  consistent  set  of  length  units  is  necessary  - 
inches  or  feet. 

PROGRAM  CAPABILITIES, 

Analyses  can  be  performed  for  surfaces  with  or  without  cambering,  and/or 
with  of  without  a  control  surface.  Analyses  can  be  performed  for  1/k^  ■  0  to 
09 .  ±/k  =  0  is  equivalent  to  zero  forward  velocity,  and  yields  the  aerodynam¬ 

ics  associated  with  a  ground  vibration  test;  1/k  -«  «  is  equivalent  to'  the 
steady  state  fiow>  case,  m  =  0.  The  surface  may  be  divided  into  as  many  as  25 
strips.  The  number  of  reduced  velocities  used  in  any  one  analysis  (one  input 
deck)  must  be  50.  If  it  is  desired  to  compute  the  matrix  of  steady  AIC£ 

%  a  negative  value  of  1/k  should  be  supplied  to  the  program  (S  and  c 
must' also  be  provided).  r 

DATA > DECK  SETUP 


1.  Title  Card  1 

2.  Title  Card  2 

3.  NCAM,  ISZ,  JSZ,  N0PUNJ 

4.  cosA,b  ,s,S,c 

5.  Ay,  b,  5--,  for  each  strip 

6.  1/.  series 

k 

r 

INPUT  DATA  DESCRIPTION 

1.  &  2.  Title  Card  1  and  2  may  contain  any  characters  desired  in  Column  2 
through  72.  Column  1  should  be  blank;  Characters  on  these  two  cards 
appear  at  the  top  of.  the  fiyst  page  of  printed  output. 

3.  Control  card  (Format  1814) 


Column 

/:  1-4 

•  5-8  . 

9-16 

,17-20 

Name 

/NCAM 

ISZ 

JSZ 

N0PUNJ  ' 

Item  ; 

(1)  '5 

(2) 

(3) 

(4) 

NCAM  =,0j  Camber  not  considered  in  analysis  (rigid  chord). 

=*'  1}  Camber  included  in  analysis. 

ISZ  ~  Number  of  strips  25 
JSZ  •=  Number  of  reduced  velocities  <50 
N0PUNJ  =  Or*  1_  &  |  Ch]  Punched  out 
^r 

l~No  punched  output 


(13) 


4.  Data  Card  (Format  6E12,8) 


Column 

1-12 

13-24 

25-36 

37-48 

49-66 

Name 

C0SLMD  - 

BR 

S 

CAPS'  ; 

CBAR  _ 

Item 

(1)  ~ 

(2) 

(3) 

(4) 

(5) 

CjSSLMB  *  cos  A 

BR  =  b 

r 

S  *  s 

CAPS  =  S 

CBAR  *  c 


*  cosine  value  of  the  quarter  chord  sweep  angle 

*  reference  semirchord 
=  semi-span 

*  surface  area  of  wing  (required  only  for  steady-state  analysis) 

*  mean  aerodynamic  chord  of  wing  (required  only  for  steady- 
state  analysis) 


5.  Data  Card  (Format  6E12.8)  Repeat  for  each  strip. 


Column 

1-12 

13r23 

25-36 

37-48, 

.  49-60: 

- - - 

81-72 

Name 

DELTAY 

B 

■/;  V  . 

:  Z2 

Z3 

Item 

(1) 

_<2>  . 

(3) 

(4) 

(5) 

DELTAY  =  Ay^  =  Strip  width  of  strip  "i". 


B  =  =  Local  semichord  of  strip  :'i". 

Zl  =  =  Fraction  of  chord  for  location  of  forward  control  point  on 

i  strip  :,i".  When  NCAM  =  0,  must  equal  .25. 


Z2  ■  ?2  *  Fraction  of  chord  for  location  second  control  point  on  strip 

i  "i",  Is  negative  for  control  surface  on  strip.  When, 

NCAM  =  6,  ?2  oust  be  the  percent  chord  that  corresponds  to 

the  control  surface  hinge  line. 


Z3  =  S3  =  Fraction  of  chord  to,  third  control  point  on  strip  ’’i".  When 
i  NCAM  =  0;  Ernst  equal  zero.  When  NCAM  =  1,  must  be  the 

percent  chord  that  corresponds  to  the  control  surface  hinge 
line. 


NOTE:  When  NCAM  =  0  and  ^  is  negative,  is  located  internally  in  'the  program 
and  is  placed  at  the  trailing  edge  (g^  -  1.0). 


The  distance  between  the  forward  and  middle  control  points  (d)  and  the 
control  surface  local  chord  (c  )  are  computed  internally  in  the  program  and 
are  printed  along  with  the  strip'  data  in  the  program  output.  When  NCAM  =  1 
and  ^2  Is  negative,  is  located  internally  in  the  program  and  is  placed  at 

the  trailing  edge  (t^  =  1.0).  The  control  surface  local  chord  (cfl)  also  is 
determined  internally  in  the  program,  and  can  be  found  in  the  printed  output 
with  the  geometric  data  for  the  strip. 
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6'.  Data  Card  ~  1/k  series  (Format  'SSI 2.8) 

Column 
Name 

Item 

1/k  *  1/k^  *  Reduced  velocity  V/b^w;  continue  oh  successive  cards  until 
i  •  JSZ  (6  values  per  card); 


1-12 

13-24 

:25-36 

37-48 

49-60 

,  61-72 

r/k^ 

l/k2,  ' 

:  t  -  ' 

1/ki 

i 

•  «  • 

1/kJSZ 

; 

o 


I 

l 

I 

1 

I 


> 
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EXAMPLE  PR013LEM 


As  an  example  problem,  the  subsonic  AIC's  are  calculated  for 
the  high  aspect  ratio  svept-back  wing  shown  below.  The  wing  is  analyzed 
for  the  rigid  chord  case  (no  camber)  and  for  the  flexible  chord  case 
(parabolic  cambering):.  The  analysis  is  performed  ior  the  reduced  frequencies 
(1/k)  of  0,5. 0,8.0,  and  -16.0.  A  1/k  -  0.0  calculates  the  aerodynamics 
associated  with  a  ground  vibration  test.  A  negative  1/k  calculates  aero¬ 
dynamics  associated  with  steady  state  flight. 


Note:  Any,  number  may  be  used  for  the  steady-state  case  as  long  as  it 
is  negative. 


3 


niCGRAM  INPUT  DATA 
NO  CAMBER  CASE 


■i 


cosA 

b 

r 

s 

S 

c 


*=  .75 

=  6.0  ft 

20.0  it 
=  200  ft2 

=  15.0  ft 


1/  =  5.0,  8.0,  0.0, 

k  h  257  chordline  for  the  no  camber  case.  t2  must  be  on 

H  must  aJ"ayS  t  '  surf3ce  is  present!  *en  bo  =?»trol  surface  Is 

hinge  line  when  a  control  surface  .  . 

present,  t2  «*  »«  W  P°Slti°”' 


NOTE: 


1 
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PROGRAM  INPUT  DATA 


STRIP  NO. 

AY 

Kft) 

'l 

52 

1 

3.8 

7.5 

.3 

.55 

.803 

2 

3.6 

6.8> 

.3 

•55 

.779 

3 

3.4 

6.2 

.3 

-.55 

.734. 

4 

3.2 

5.5 

.3 

-.55 

.723 

5 

3.0 

5.0 

.3 

.55- 

.700 

cos  A 

11 

• 

b 

=  6.0  ft 

r 

s 

=  20.0  ft 

S 

=  200  ft2 

c 

=  15.0  ft 

1/k 

*  5.0,  8.0, 

0.0,  -16.0 

£  1* 2*  an^  S3  ^  ba  any  .arbitrary  position  when  no  control  surface  is  present^ 
when  a  control  surface  is  present,  ? ^  and  ^  raay  be  arbitrarily  located  and  ? ^ 
must  be  located  at  the  hinge  line.  In  both  cases,  however,  ?2’  anc*  ^3 
should  be  distributed  across  the  chord  so  that  the  chamber  can  be  properly 
defined,  e.g.,  =  .20,  ? ^  =  «50,  and  ? ^  =-80. 

NOTE'  Negative  ^  for  strips  3  and  4  indicates  a  control  surface  on  these  strips. 


(18) 


uj  u:  u>  ui  as 

—  IA  ^  W  ro  cr 

w  O'  &■  o  d>  o- 

Q.  Wr<  O  WW 

N  O  J\  ▼ 


3  3  O  3  3 
c?  ooo 


X 

U  UML  LL  X 

1L 

— 

d  '3  03  d 

3 

05 

—  ■**-  CM  ’*5  A. 

** 

3-3  3  <3  d 

3 

-Odd 

cm 

-5  ^  •«■  fO  3 

3 

X 

' 

tr^ 

C^'B’W  c. 

3 

— J 

a.  ii lXj  o  u- 

WA- /*■*>.'  tv. 

•W 

.‘X 

cvd  e>  o  -o 

•4  >4  <4  >4  4 

«4 

— 

03 

o  o  =  e:  o 

1 

3 

,i,*5 

05- 

oo  3  o  e 

u. 

X 

M 

-Td  rd^* 

X 

X 

•<. 

"*»  -G'SV.CM  -r-l. 

X 

.X 

vw 

•♦4  4-4  • 

•  O  3  =  •  C  d 

"< 

-AT 

o 

— :  ^  —  —  — — 

—  •— •— 

Ci 

o 

X 

C5 

•IT 

H 

>#■* 

J.  -Li-  1U  *li  u. 

X 

-X 

*— 

g 

It  II  II  M  If 

>*— 

o  e  e  e-.e. 

■X 

-— 

X 

w 

3  3  3  3  3  ■ 

>x 

X 

-X  - 

X  -<C  X 

r-o 

3  3X33 

— 

•>- 

7 

JJ-Z 

r.  <  x.  < 

N‘ 

x  it  x  r  r 

1/5 

X 

>u  x 

X  o.o.  it  ,x 

CM  CM  CM  CM  CM 

r 

—  3 

x  x  ur-r  - 

■«  •»  >4  4  4 

■4 

X 

tr 

a:  c 

<^L?  1/  C 

d  c  o  trr  e 

3 

o 

— 

x 

-1  1  —  X 

2. 

.1/)  a: 

—  XX 

— 

__ 

X  d.  -*.  Q 

♦— 

o 

*— 

-T  xr 

-Z  IL  CC  -U 

X 

.X 

'(J 

—  o/  a: 

T- 

X 

•— 

«A  X 

0-3  0  3-3 

3 

■a. 

-i/; 

U 

o  a.  v> 

a 

u 

U 

c  U 

#■» 

>u  ll  a  u  a 

LL 

a: 

X 

•x 

.a: 

— • 

0  3  3  3  O' 

d 

X 

c. 

ll 

■w 

3  3  3'3  3 

c 

>• 

•X 

X 

X 

3  3  3  3  3 

— 

o. 

try 

LL 

IT.  X-  CM  IA  C 

3 

c. 

— 

_ 

N  £  tx  X> 

d 

LL 

t r. 

LX 

IL 

{4 

y 

>- 

z; 

V- 

—  ■—  —  — 

— 

_  ■ 

.X 

■•< 

X 

>a 

2 

. _ 

«— 

-1 

a 

f— 

•v 

— 

cr. 

O' 

*c 

w 

-A 

—  T—  —  -«-4 

r- 

-5JJ 

o 

'*— • 

3  3  C  3  C. 

3 

— 

.w 

•y 

>W 

<< 

< 

L  L  ll  U.  Li 

u 

> 

I 

> 

d 

2L 

> 

C  3  C  C  C. 

a. 

X 

X 

< 

_  L  .,  „•  ^ 

o 

a 

-*5 

/X 

H- 

(T  CT  M  O 

d 

U 

X 

_■ 

n  r.  n  n  n 

iA 

«x 

.‘X 

X 

•  4  '4  ■•  • 

>4 

•C 

*< 

X 

X  3  3  3  X 

d 

•*-*  *H  O  rtO.O 

GOO  OOD 

LU  LU  LU  Ui  Ui  Ui 

▼  n  3)  co.-o  o 

n  <\J  h  co 

o  >c  o  m  eg  o 

co  oo  v\ ▼  a> 

rt’ON  <u  ■<-  rv 

fO  CM  ▼  O  *3  *0 

C3N  fO  3>  CM 

•CHjy  n  CO  CM 

•  •  •  «  • 

CD  O  O  C  OO 

II  If 

N  H  fl  rl  H  fi 

o  —  g  a'£?*=> 

u:  m  n  i  uj  Lb 

X  .>  CM  in.  **  fv 

oca  ▼  o  < 

23  UV  'O  O'  *C  CM 

CM  O.  !*?  G  **  fv 

rv  cmg  *-<  uv  re 

O  O  tr»  H  T  03 

o  o  »o  .-v  re  cm 

if*  *—  air  r< 

*  •  •  •  •  • 

G  W  G  .=  G  G 

III  III 


CM  r-> 

O  r* 

n  CD 

C 

HG  rH 

r-i  r*i 

=2  23 

=5  — 

P  ° 

~ 

GOG 

o  o 

1 

a;  —■ 

-b  a. 

a?  u-  ui 

Lb  Ui  uj 

LU  lb 

o  X 

X  *-* 

o  o- 

V  CM  X 

e  CM 

P  P 

o  o* 

m  cm 

»e 

O'  .-V 

ro  rH 

r-  •j' 

CP  re 

re  *c 

fv 

■OW'D 

w  irv 

f* 

** 

CM  *e 

> 

O  i 

P  X 

e 

jk  » 

CM 

e  g 

o 

!f\  K 

cm  rv 

'O  *T- 

r-» 

*c  x 

o  X 

>- 

CM  P 

>- 

4*0  <r 

D- 

C' 

M" 

>- 

rv  *o 

>- 

o  ^ 

Lb 

II 

X 

C  CM 

C. 

w 

X 

O'  vC 

CM 

X 

rtsc 

a 

f*3  X 

.=> 

CM  CM 

▼  CM 

▼  rt  O' 

H  H 

CD 

cr. 

ev¬ 

•  • 

CM 

«  •« 

n 

•  « 

• 

ie 

•  •  • 

CM 

j  • 

» 

X 

O  S 

w  X 

P  G 

O 

o:p  o 

O  P 

o' 

— 

il 

1 

II 

1 

ii 

1  1 

n 

1  1 

II 

1 

ir\ 

h- 

Lb 

CM  O 

U. 

CM  G 

iL 

g  ri 

rH 

.u. 

Or  fi- 

LU 

CM  o 

• 

V> 

fM 

X  G 

JV 

2  — 

N 

G  o 

o 

.VI 

CS  O 

M 

=  o 

<= 

»- 

•— • 

►— 

u. 

u: 

Lb  lb 

tr. 

lb  -b 

w 

Lb 

tn 

a.  u  a 

y; 

Lb  ib* 

ii 

o 

®  cc 

▼  re 

X  X 

CV. 

IO  CD  CV 

X  <= 

C  0* 

— * 

*-«  sn 

— * 

*3  fv 

X 

*** 

e  x'  *c 

<-» 

O'  CM 

X 

r- 

rv  o 

CM 

G  X 

y 

C  P 

CM 

V 

f*ON 

IT 

CM  X 

X 

r-t  C 

o  p 

»\  X 

T 

XI  fv  .0 

M-  CM 

w 

w* 

IT  *■ 

■W 

'C  ^ 

w 

O'  -  o 

>C 

w 

e  cv  ^ 

w 

M*  >0 

X  » 

O'  .s. 

i  .v 

fl  r*  if 

X 

y  rO 

N 

to 

V  r* 

u 

V  CC 

c 

X  X 

rv 

to 

If  CV  ^ 

to 

IfV  r- 

r* 

X 

ri  CM 

CM  O* 

CM*  O'  r-* 

n  n 

•  • 

•  • 

•  • 

• 

•  •  • 

•  • 

1 

7  ~ 

-  G 
1 

w 

n 

i) 

(i- 

1 

^  <“ 

•-4 

—  r- 

T-»  f-<  «“< 

»-f  r- 

G  G 

S  G 

2  c  c 

e  g  g 

O  G 

U,  Lb 

a  u. 

a  u.  u. 

Lb  Lb  lb 

lb  lb 

CM  T 

e  - 

a>i  v  — 

v'  -*v  V 

=  CM* 

re  o 

re  O' 

O  O  fv 

CM  re  G 

o  n 

CM  X 

c  *e 

X  p  x. 

O'  r-i  G 

v*  X 

re 

X  <r 

rv  m  «~ 

re  *c  g 

re  g 

x  > 

x  iv 

"V  Pf 

rv  w  cm 

fv  x 

rv.  aC’ 

TT 

o  xi 

r-*  O  r- 

CM  W 

X  CM 

X  *- 

N  .\  «C 

'C  r  CV 

vC  X 

re  cm 

re  cm 

r  r  N 

r4  »-  e 

re 

•  • 

•  • 

•  •  • 

•  •  • 

•  • 

G  G 

1 

7  = 

G  G  G> 

1 

G  G  G 

1 

G  C 

1 

C' 

A-  _ 

^  —  G 

—  *2 

CM  r' 

G  G 

re  g 

G  G  G 

G  C 

lb  LU 

Lb  LL 

U.  LU  lb* 

Lb  Ui  UJ 

UJ  lb* 

X  CM 

M  CM 

O  30  r-1 

CM  -X?  O' 

CM*  CM 

O'  X* 

e  x 

M  L*V  a 

O'  V  G 

CM  M" 

•e  iv 

r  e 

C  C  X 

=  NO  O- 

M-  O 

X  X 

G  X* 

<.  re  o 

C  H  C 

O  Xv 

O'  rv 

X  M- 

r  ■m  o 

rv  o  X 

O'  CM 

f—  r-1 

re  o 

O  IV  IV 

e  t-p  x 

0  rv 

M  P 

e  <r 

rv  p-  e 

iv  X' 

e  cm 

r~  r“ 

f~  ,*v 

o  r.  ,e 

o  re  re 

r-  CM 

•  • 

•  • 

•  •  • 

•  •  • 

•  • 

r  c 

!"  *C 

•  t  — 

.  ••  c 

o  c 

II  II 


(DUJUJ 

UJ  LU  UJ 

£  O' 

to  O  A 

tir 

o  irv 

fo  o 

r-l  eo 

aik 

O  U>  •«• 

to  A  to 

Ok  so 

QC 

r-  o-  n 

A  A  CO 

ai 

• 

lA  IV  ♦ 

O  A  Ak 

*** 

cr. 

o  n 

O  V*  CM 

r 

•  •  • 

•  •  • 

< 

O  O  G 

o  o  o 

1  1 

t  1 

— 

- 

oj  o;  ih 

NCJH 

my> 

a 

G  C  S3 

coo 

H 

3 

JUJU! 

UJ  UJ  UJ 

h 

p* 

kf\  X  X 

r*4  X>  rkk 

A 

•— 

N’lO  IA 

t-i  X  o 

JC 

0‘v  in 

JO  to  CM 

oj  «o  o; 

O  CM  k© 

3- 

.o  a  v 

rtiA  C 

+•0 

X 

OON 

■O  A  O 

Jj 

w> 

Oy  *  *\ 

<e  :o  « 

if 

UJ 

CM  r-  *0 

CM  T-»  to 

X 

X 

•  •-  « 

•  *  • 

r— 

o  c.  o 

O  C  Ci 

1  f  1 

1  1  1. 

a 

— * 

'?• 

a; 

f.  t-i 

r* 

CM,  G  O 

H  rt  H 

i’ : 

— 

S  <3 

■S3  O 

OSS 

C  G  O 

o  o 

A 

Ot 

- 

1  1 

X*  OJ 

OJ  JLi 

X  Uj  LL. 

UJ  til  Ui 

Ui  LU 

a. 

N-  r%k 

i*. 

^  ▼  CM 

008 

CM  C 

-j 

t-l  A 

•m  ,-vj 

V  ^ 

O  T  O' 

O  CM 

pm 

X 

r*k  ** 

**0© 

G  G  + 

cc  «o 

ft  X 

*? 

— 

ft 

A 

Cv' 

3  rs. 

»\ 

o 

O  T 

>c 

Ok 

CM 

»-l  X 

Ci 

w 

G 

>C  o 

T-  A 

ri  « 

CM  to  ▼ 

tH  CO 

w 

>- 

=3  A 

>- 

cm  to 

>* 

M3 

>■ 

CM  CVJ 

>- 

CM  G 

u: 

UJ 

II 

t-  X 

c 

<.  T 

c 

G  V  Ok 

X 

^  CM  ^ 

X 

f>  o 

5 

C* 

w 

tn  .o 

*■' 

n  .o 

H 

J'NO 

VT'CM 

< 

X 

ci 

C/J 

C\J 

•  * 

cm 

•  • 

to 

•  •  "• 

to 

•  e  • 

CM 

•  • 

X 

s 

0-. 

c  cr 

O  =J 

GOO 

o  o  o 

s  o 

o 

S3 

•-* 

II 

1 

fl 

i 

II 

1  1 

II 

1  1' 

II 

1 

It 

o 

S 

X 

^"t“ 

z 

X 

K 

u. 

cm  c 

X 

Cyo 

X 

r»  W  n 

LL 

t-i  CM  tM 

UJ 

CM  o 

w 

“ 

• 

»* 

N» 

—  SS. 

fsi 

w  w 

Os. 

SOS 

■N 

GOO 

Kl 

o  o 

c 

M 

*— 

*- 

k— 

•— • 

a 

>- 

a 

In 

X  X 

tn 

X-  X 

cn 

\X>  LLi  vJJ, 

a*  uj  uj 

t/> 

uj  uj 

X 

II 

w 

£  X 

in  to 

X  0--0 

CM  O  CC 

rx  o 

f 

—» 

OO 

*>4 

X 

O  NO  «- 

*—> 

•O  CM 

1 

(A 

f 

X 

t-» 

X  o 

c\ 

cr  cc 

r' 

to  *r  to 

V 

IP  O  c 

A 

X'CC 

1 

t— 

X 

Uj 

•x  o 

V  0“ 

O  CM  A 

t-,X  t- 

•O  CM 

m 

2 

tf 

a 

Cs?  V 

W 

Pv  ▼ 

k» 

«CN 

w 

O  AO 

w 

CM  O 

A. 

•1. 

— i  *r 

X 

.o 

X 

c  to  ^ 

X 

■c  ,-o 

> 

X 

ti 

fk  f 5 

Ci 

1-*  X 

L* 

X  IA  ^ 

c: 

CM  to  CM 

u 

CM  ▼•< 

w 

r- 

X 

C  CM 

V  r< 

CM  OJ  V 

CM  CM  V 

t-i 

m 

■— * 

•  • 

•  • 

•  »  t 

i  ♦  * 

•  • 

g 

u. 

3  X 

X  X 

GOO 

o  o  o 

O  S3 

fl 

L. 

1 

1 

1 

M 

kV 

C 

o 

«— 

H 

«k  «-( 

rV  •— i 

H  H 

c  o 

o  c 

GOG 

GOO 

O  G 

n 

fj 

1 

1 

9 

w 

x  a 

a*  x 

a  IblL 

UU.  U 

UJ  UJ 

| 

” 

a 

x  r- 

VC  X 

r  c  c* 

o.  o 

■ 

UJ 

O'  a 

X  cv 

r?  x  »-* 

f*  CM 

X,  CM 

m 

JT> 

■c  r*» 

X  c 

Ok  o  if; 

M*  A  A 

*o  X 

_ 

o 

c  cv 

3  <N 

fO  «f\s  IO 

CM  X 

U 

*c  o 

OJ  x\ 

a  r  o 

O  Cv  *o 

A  X 

X 

C  A 

to  ?0 

O*  CM  CV 

a^s 

t-k  e 

■ 

— 

, 

o  c 

Cf  c 

r-  r-  X 

Tv^X  X 

<  O 

I 

a  to 

O'  *o 

^  CM  X 

ICrt  ^ 

X  CM 

I 

c 

•  • 

•  • 

•  •  • 

•  ♦  • 

•  « 

■9 

•3  S3 

CS  G 

o  cr  c 

o  o  c 

o  c 

3. 

1 

1 

1 

1 

1 

<* 

' 

^  r. 

O  — ' 

o  —  r 

0  C 

c\  r- 

I 

> 

x  e 

C  X 

o  o  cr 

COG 

cr  o 

I 

c 

x  a 

X  lL 

UJ  u.  u. 

UJ  UJ  UJ 

UJ  U* 

1 

CL. 

hs  cm 

X  OJ 

CM  M  O 

-0  to  *> 

O  CM 

11J 

r*k  ir 

O.  X 

ro  oN  cc 

H  V  CO 

v  vr 

=  fs 

1-  c 

•A  X  X 

O  O’  • c 

o  c 

CM  X 

if\ 

nc  r>.  c*> 

O'  M*  A 

vr  a 

* 

so 

to  rr 

A  O  f*k 

Ok-vr  t- 

t-t  CM 

I 

X  T-> 

o.  a> 

V  S  UN 

lf  V  c 

3 

o  =. 

o  « 

<  x.  — . 

A  O  O 

t-i  CM 

1 

tO  f 

rr  fk 

Ok  c  o 

cv  o>  cc 

vr  cm 

»  • 

•  • 

•  •  ? 

1  •  • 

•  • 

■  «* 

-  5T 

•;.  s  - 

«r  -S  — 

C*  w 

t  i  r  7 


i 


[ 

r 


CVJH 

CUv-f.CSJ 

75  C  O 

O  O  C3 

1  (  1 

1  C  1 

Ui  IU  UJ 

:>  s>  in. 

N  CC  O 

>C  'C  o 

o  ro  o* 

CNO> 

o  «c 

cnfi 

croo 

ri'C  O 

▼  in  cu 

r4  W  f 

cv  cv 

CO  'S?  O 

CO  TfO 

CC  *-  *-l 

>0  n  CD 

f. 


•  «  ••  •  •  •  •••  •  • 

©  ©  ©  ©  e  ©  ©  coo  cc 


r  ^  c-  f  :\  -  — 

©  e  ©'©  ©  c  c  c 

I  t  I  II 


C 

U.  LL 

IL  U. 

ILl  LL  UJ 

U.  LL  LL 

U/LL 

a 

ri  cv 

rs  cv 

W  O  X 

cv:  cv  :o 

X  CV 

a*  ^ 

co  in 

©  © 

llUO  £► 

fs.  o 

m  v 

■ 

©  fs. 

O'  <V  O 

J  HX 

**  'C 

1 

cc  cc 

<■  lf» 

x  e  x 

xn^ 

fv  lf» 

w 

v  r* 

ro  v 

iT  =  c* 

N  V:  fO 

in  cv 

fH  fH 

cv.  o 

*0  cr  t-i 

O'  O'  CV 

ro 

X  c 

•-*  v 

CiT  t 

r*  O  CC 

'C  cv 

#  CV  r* 

cv.  r^ 

r*tf  CC 

'C 

co  cv 

m 

•  • 

•  • 

•  •  • 

•  •  • 

•  • 

1 

©  <r 

cr  c 

IT  C  7 

C22) 


u 


*  n 

*■  CM 


cr 

x 

< 

15 


>* 

cr 


I 

fc  I 


o 

f-  c  o 

fH  O  O 

u: 

* 

coo 

0-0  0 

-,‘X 

'r— 

Ur  U.  Ur 

Ui  Ur  lb 

X  CZ  UN 

n.O*  n* 

Sr 

n  CV  r>. 

UN  O  x 

** 

IV  M)  XV 

xn  cv 

X 

UN  O'N 

m  rs..^ 

w 

X*  O  >* 

m  to  o* 

-W 

o 

M3  *+,r* 

T*  >C  ** 

^  x  <n 

CO  H  K 

a. 

II 

r-  M3  *-4 

H  <  n 

-j 

•  '•  • 

•  *  • 

a 

*< 

o  o  c 

o  c  o 

* 

CD 

vO 

CM 

IV 

"5 

1  1  1 

»0 

1  1  ’ 

CV 

w 

IU 

_ 

uO 

r 

'>• 

>- 

V* 

>- 

a 

c 

il 

tt 

a 

CC 

a 

c 

x 

ri 

o  «-<  o 

o  HO 

a. 

UJ 

UV 

CM 

O 

cv 

e 

fO 

o  o  c 

to. 

coo 

cv 

c 

-X 

WJ 

« 

a 

c 

«< 

Ui 

11 

u: 

ti 

u: 

II 

Ur  UI  UJ 

II 

Ui  UJ  UJ 

II 

UJ 

t: 

o 

X 

-r 

IV 

i-* 

X  C*  *H 

cv  o  o 

▼ 

z 

X. 

*- 

Ur 

r*. 

Ur 

a 

Ur 

cv  ^  in 

Ui 

X  IO  « 

IL 

~ 

>- 

w 

x 

.■M 

cv 

Nl 

r 

KJ 

V  V  = 

rsJ 

CO* 

rv 

X 

>c 

—• 

T— 

i— 

it  cr. 

nio  -c 

->• 

< 

u 

*/> 

X 

to 

-T* 

O’ 

*-l  CM 

in 

N  N  ▼ 

to 

-c 

CL 

u> 

UI 

o 

CM 

c. 

o*  m  o 

CNO 

in 

4 

CD 

V 

z? 

T 

r~» 

CV  r-  M 

-  ri 

•n  c  c 

« 

cr. 

cn 

or 

»• 

cv 

rv 

CV 

ro 

tT'*-*  CV 

'T 

T  CV 

’X* 

cv 

►— 

< 

Ur 

•  • 

•  « 

•  -  •  * 

‘J  • 

• 

z 

-J 

a 

V* 

X  o 

w 

CO  c 

o  o  c 

w> 

boo 

w 

o 

» 

x. 

J. 

1  ' 

X 

1 

X 

X 

X 

1 

*“ 

ii 

a 

u 

o 

c 

u 

w 

2' 

» 

- 

U- 

'oc 

—< 

n 

C  r» 

— 

•*  e  r* 

r- 

u 

c 

c 

COO 

c  c  c 

c 

X 

1 

"  1 

c 

v. 

Ui 

UI 

Ur  Ur  UI 

UI  Ur  Ur 

UI 

15 

CV 

un  cv  cv 

iT  ■=> 

NT 

y, 

O 


2 


O 

.X 

lb 

<• 


CM 

'C 

X 

CM 

M* 

CM 

•  • 
o  c 


(* 

c 


NT 

cv 


rv  n*  »e 
•c  x  r- 
^  w  v 

«C  — <  ’J' 
.■‘JOW 
’*C  NX 


o  c  c 
<1  I 


▼  O'  w 
CM  tIN  O 
X  CV 
V"  r-”  CV 
'O  ’H  tf' 
M:  r^ 
f-T  r. 


CN 

T*4 

>T 

UN 

IO 

CM 

•  • 

c  c 


1 


f.  <-/ 

'< 


!l 


I 


r»'ci 

Lit 


A 


b\ 

I 


1 

I 

I 

! 


- 

H»4H  3 

•4HH3 

°?°° 

o  o  o  o 

#- 

UJ  UJ  UJ  UJ 

UJ  UJ  Ul  UJ 

^  >  in  S 

O  CV  O  A 

- 

,  ' 

oO'Nn 

«  A  CV  A 

•V'A  -o  o 

H  t’jS  O 

< 

ONtSNC 

Oho#* 

©NON. 

«  UNT  A. 

” 

»NlftT 

VN»«' 

w  evjo-  o. 

O  A  «  CV 

NNNCM 

o  -o  CV  CV 

•  '  •  •  • 

•  •  •  # 

o  o  o  o 

o  c*  o  o 

z 

I  1 

1  I 

" 

:v  cv  cv  mi 

CV  TV  cv  ^ 

- 

o  o  o  o 

o  o.o  c 

’ 

US  UJ  UJ  UJ 

UJ  UJ  MJ  UJ 

$ 

o  cv  t>  rv 

O  CO  CV  A 

o  oC)  0> 

ro.o  x-  x 

A  ^N  ■© 

TMION 

r> 

r-  u>  x-  n 

ICN  <V  A 

O  A  O  o 

■o  x  y>  "O 

^•TIO'/V 

ffi  cv  X*  «? 

o-tHn 

cv  cv  -*  cv 

H  rJ  H 

H  H 

- 

•  •  •  • 

•  •  •  • 

o  o  o  a 

o  o  o  o 

_ 

- 

1  1  1 

1  1  1 

a: 

* 

_ 

jj 

x 

CV  CV  *- 

M  O.  »< 

H  =  r*o 

CV  CV  tv 

X 

o  o  o 

o  o  o 

0  0=0 

O  C  O  CJ 

O  O  O 

-* 

cu 

JJ  UJ  JJ 

US  111  UJ 

UJ  UJ  UJ  UJ 

Oijsx  Ji 

If  UJ  Ui 

tr  un  a 

C(MH 

CV  CPT 

NCV  o 

- 

ir  a 

X 

X  A  VN 

-=  X-  3 

x  -<  u\  r«- 

JO  A  A:*c 

A  CV  A 

*— 

C9M 

r  un  to 

C  VN«C 

PO  =  ©  A 

X  -©  x- 

•“ 

•t-on 

-O  ,=  X 

r»  3X  X 

un  ^  cr  rv 

VN  A  CV 

;x 

O-  0*^0 

A  'C  fN. 

A  >6  f  I 

A  fO>X 

X 

•C  A  X 

o  ;n  r  *■ 

X  X'= 

—  A  p: 

> 

sr>  cv  x 

n  x\<c 

A  A  10  C 

X  »- 

A  A~~- 

Hr!  rt 

rr  —  cv 

— i  rv 

"3  ««*  ^  3, 

CV  CV  — 

* 

s« 

•  • 

•  *  • 

*  «  •  • 

•  ».  • 

a  •  ♦ 

a* 

CC  3 

css 

C  3  1  C 

O  ■=  ='  = 

=  o  o 

. 

x 

| 

1 

1 

1  1 

.1 

t— 

CV  PJ 

cv  cv 

N  N  M  H 

CV  CV  CV.'«-» 

CV  CV  CV 

0. 

o  o  s 

r  c  3 

' 

=  «=,  0.0 

X  X  =  = 

o  =  ■= 

a 

-U.  U!  If 

UJ  JJ  JJ 

- 

UJ  UJ  U  uil 

UJ  UJ  il  UJ 

JJ  if  UJ 

*- 

't‘w  UN 

cco 

•Ch'ON 

.ox  x  e 

CV  X  A 

vi 

▼  -3  X 

—  a  — 

o  JN  3  -  *0 

▼  to-M  rv 

H  S  *1 

CN  ■©  CV 

a  r:  a 

X'=  O  O 

•C  X  CV  A 

—  CV  ifN 

ai 

•x  37  CV 

■o  =  x- 

x-  cv  *r  rv 

CV  «*■*  -\t  -1 

*r  —  cv 

ClfiN 

«.uv  a 

flO  CV  ffi  X 

c  s  Cv  — 

A  •—  “3 

X 

X  "3.3 

jo  a  o 

«  -•  O  X* 

A  x  35 

•XN  *o  lO 

•— 

r"  (V  <0 

IC>  NT  *-* 

rt  Ti.TN 

CV  CV  *- 

SOX  X 

n 

cv-o  — 

re  •■«  cv 

H  JO  «  r< 

H.'iH.N 

Pro, 

CO 

•  • 

*  •  • 

•  •  *  • 

•  •  •<  * 

••  * 

,1* 

TN 

.o 

X  o  X 

**5 

3  3  = 

X- 

=  03'= 

•c 

0  0  =  0 

'*■* 

=  =  = 

a: 

S 

r  i  i 

1  1 

1  1 

-  1  -  1 

i  i 

CL 

>- 

>- 

•>- 

y- 

>-. 

X 

Ut 

u 

ffi 

CT 

cc 

z 

e 

c 

O- 

=5 

Aicwh 

:v  cv  -» 

• 

rv'W  — 

CV  CV'—  = 

£v  sv  cv 

■J 

c. 

CO 

O  o  3 

7? 

=  03 

X* 

o  =  =  c 

-T 

=  o  =  c 

yj 

z 

3 

X 

— 

C 

II 

LL  UJ  UJ 

11 

u;  ll  uj 

II 

u:  UJ  U?  UJ 

M 

UJ  u.  UJ  u. 

II 

UJ  1*.  UJ 

= 

, 

o  >-"= 

-3  =  e  a 

A-  .O  A  » 

O  A  O 

>- 

<IT 

►*- 

If 

NCH 

IL 

COO 

;ix 

<N«r/ 

u: 

o  cv  CtA 

UJ 

A  X  O 

r 

• 

vO 

VI 

c«c 

A 

V*N  **3  A 

IM 

O  -3  x  =, 

A 

rv  \r  o  a 

A 

'=•=  o 

3 

~- 

Ht- w 

mm 

*«*- 

mm 

T  C  T  r< 

mrn 

X*-  CN  tO 

*  mm 

xic.tr 

ir. 

CO 

a  «o  o 

V3 

'=  P  X 

Vi 

>  X  X  X 

Vi 

—  A  C  'P 

cv  c  c. 

— 

W 

X, 

a  ip 

X  UN  «o 

A  ’^n  r-  f- 

If  s«x 

AAA 

, 

z 

— 1 

•JV  ON 

*** 

A  ^  =•. 

C  UN  cv  o 

A  —  m.  — 

?  r  p 

tr 

cr 

X  —  10 

cy 

CN  *-  (T 

r- 

•**  o  ^  ^ 

V 

»-  —  *0“A 

6* 

*7*0  — 

— 

X 

X 

•  ’•  '• 

•  •  • 

•  •  •  • 

*  •  •  • 

•  •  • 

o 

w 

CD 

O  C3  O 

«— 

=  c  o 

w 

=  0=0 

w 

=  0  =  0 

w 

0  =  0 

— 

JL 

X 

1  1 

x 

;ti  *  ;r 

.i 

I  II 

X 

1  t  1 

1.  1 

u 

V 

X 

c 

c 

u 

’ll 

u' 

u. 

r* 

X 

tv  »\  cv. 

—  CV  CV. 

rv  rc  — 

■  N  ,N  .\  - 

cv  -- 

.11 

coo 

'3  X  C 

e  sc  c 

C*  =  c  =  - 

A  C  = 

tS 

iuia- 

a  a  l 

u  u  a  a 

u  u  u  u 

U  U  u. 

-=-  x  ▼ 

-  <■ .% 

<0  O  JN  = 

o  r  r.  a 

—  A  X 

u. 

nn  irv 

X  c'x 

X*  X  X  A 

Atn3 

X  c  — 

'-5, 

K.  un 

cv  *o  = 

.■V  A  o  JN 

73  O  X  a 

z- 

nr  un  « 

XXX- 

CCNOO 

’X  A  X  -CN 

C  «*.  CV 

N  x  n 

cv. 

CN  A  .V  -IN 

•C  —  =  fO 

—  — 

p*  x  ir 

<  tr.  tr. 

tr  <■  <  r- 

CV  IP  A  X 

C  <IC 

-J 

•C  K  CV- 

fv  >C  V- 

•C'C?  — 

—  TAX 

u. 

—  IP  r- 

c  ▼ 

W  C  »•  C 

■C  N  *-  <C 

r  .ir  -c 

X 

•  >•  • 

•  •  • 

♦  •  •  # 

♦  •  C* 

*  i<’  ♦ 

c  e  c 
* 

w  c,c 

c  c  v  c. 

croc 

1 

1  1 

l  I 

1  N 

X 

»-  e  *- 

—•  c 

*-  X  *-  c 

^  — 

< 

3  0  3 

AC  ~ 

3  'T  — ■•=• 

c  r  r  i 

X 

1 

> 

•L  if  uj 

a  x  a 

x  *x  X  If 

a  u  t.  .f 

U‘  *1  f 

c. 

•^atc 

tr  •-  v 

x  ec  c*  a 

x  cva  r. 

XXX 

z 

v  <  p. 

'  *.  < 

.*  ^  c  — 

C  A  V  \ 

cr 

-:CN  «• 

9  *c  < 

c:  r  r  c 

X  CN.*-  *C 

«r  =  p 

L. 

0-CV.CN 

a  O  «•• 

X  CN  A 

—  J"  A  A 

A  c  .* 

N  W<C 

uv  CV.  X 

V.  A  CN  A 

9>  A  X  A 

C  CN  CO 

- 

—  CN,X 

cv  tr'  v 

•C  A  «C  C 

=  X  X  CV 

A  A  CN 

r»3  Xv 

cc  cv  cv 

A  A  X  »- 

J'  X  I  c 

•o  A  X 

•*3  *•«  SN 

to  o.  cv. 

^  cv »-  r 

Pj  *-  *-  fO 

U  »0  r- 

- 

•  «  • 

•  •  ,• 

•  •  *'N# 

•  •  •  • 

•  *  • 

3  C *  A 

t 

cr  =  cr* 

•'  « 

C  "  1  ~ 

1 

1 

CN  <N  r* 

CN  CV.V 

:(\  CN  CN  i- 

CN  CN  CN  *-  - 

CN'CN  *- 

e  o  e 

e  c  r 

o  =  =  = 

o  c  c  o' 

XXX 

a  x*  x 

a  _  j.  , . 

L.  X  X  — 

~  a  u 

Ift  A«.A 

V  CN  ^ 

-X  X  cc  A 

c  ,X  V  pi 

x  <  *- 

/“  -c  «r 

•C  o.  <, 

X/C  n  r» 

OCX 

cr  *-  C> 

O  UN  «- 

O  CV  *O  fO 

JN  A  A  A 

x  r  ^ 

=  tr.  n 

Affix 

CN  CN  X  CN 

Cvr.iT  *- 

*-  A-  X 

CN  A  *0 

x  x  a  rc 

«  %r  *-  < 

X  X  *- 

•cap. 

RTN 

xr  <  *r  c- 

tn  t-  o  tr 

>C  s  ^ 

'CCr 

K 

A  A  *C.O 

x.ir  >c  = 

in  IP 

W  X  P3 

r>  t  cv 

tO  ^n  ri  .O 

O  IT  r-  O 

cv  r,  so 

*  •  • 

•  •  • 

• 

•  «  •  • 

•  •  • 

o  c  o 

o  oc 

-=o  =  c 

o  o  o  c 

e  c  c 

1 

1 

"  ;r  f 

•.1  1 

*  »  l 

\ 


(SJO 


.(?*> 


I! 


r 

«c 


{' 


•  •  **,  .» 


o  a  ixv t,  SOON 


uiuu  turn 
».*n  e  n* 
in^u  «  • 

r  1}  OO' 

MftWO 
KV-O  N*  * 

»'r» 0 

H-  TjO  O' 
I90ft4 
•  ••>•• 


uj  uuu  uj 
«N  0*'o  O 
•  0 
>h  0  o  *r> 
7e'TT  **4 
0.HNN 
n*n  r 
in  to  e  to 

CM  «  tO  CD 
•  •  •••* 
O  ••  13 

«■  .«:■' 


i 


? 

II 

ii 


■I 

i\ 

j 


i 


H 


1 

i 

« 

■ 

I 


>» 

* 

c 

#  •  • 

*  •  • 

•  •-  •  •  • 

•  •  •  • 

«k  «r  • 

u. 

C  €=.■» 

CSC 

t»  S  3  C 

o  o  o  c, 

X 

»— 

owe: 

•».«  «= 

-,w  «•  W  O 

W  fO-H  C 

c  CC 

c. 

S>  (S  S' 

SO  J 

cccs 

c.b  =. 

«• 

1 

' 

(X 

ai  usuf 

ux  iu  ur 

UJ  UJ  X  *11 

(U  U.IU  Ul 

Till  UJ  UJ 

*— 

▼  vo  0 

cm  w  n. 

T  «  C3  «D 

▼  1ft  IO  to 

*<0  o  »o 

O 

Pv  CM  0 

-in  o  oj 

AN  HN 

OA‘itT 

0<lft 

0  x  x- 

'  «e.  *4  ♦ 

H  O  O  r« 

TCOC 

0  H  -. 

'•li 

k:  rv  «\- 

<> 

OONH 

, 

«oh>oh 

NC\' 

n*  o  -c 

T-  *0  If* 

TO  ION 

<ohn 

c  x>:< 

zz 

n  i  j\ 

■O  n  J) 

AchO  >v 

to  b  <c  n 

TTO 

- 

tOHH 

Tv0  C  « 

rtNO-O 

<  N  CM 

<r 

w 

H  N  CM 

HmM 

HtOOT 

W  -*O  CM  CO 

> 

hMO, 

*2 

cn 

•  •  • 

•  -  •  • 

4  •**  • 

i  •  • 

a 

-D 

s  s  e 

*o 

.3  S  S 

c 

“.NT 

»  c:-=>  x 

fO 

=>«==: 

X 

t 

- 

"  '  1 

-,l  • 

•  -  i 

« 

b 

c. 

>* 

>- 

>- 

■>• 

>* 

5 

X 

11 

-X 

X 

& 

X 

- 

CC 

CM 

y 

tj 

in 

to¬ 

to 

▼ 

T“ 

to 

o. 

X 

a. 

«— 

—  - 

ll'* 

u 

N 

«: 

M 

X 

i- 

U! 

u; 

tr 

UJ 

UJ 

X 

€* 

v> 

N 

N 

IN. 

ft: 

•— 

.  *— 

»L 

•— 

M 

sis 

tk 

w* 

X" 

u; 

X 

w 

ii 

w 

_ 

x 

«*■» 

a* 

a-: 

T- 

IM 

to 

T 

«* 

*— 

X 

x 

•  •  •._ 

•  •  •• 

•  •  •  • 

•  •  •  ♦, 

♦  •  » 

X 

'X 

coo 

%* 

o  o  o 

k# 

O  O  O  O 

•» 

e  e  o  c 

;c  o  o 

— 

z 

X 

dL 

X, 

x 

" 

X 

*’ 

X 

N 

.X 

C2 

C 

o 

c 

« 

•  w 

'*• 

s  s  *n 

JL'C  w 

*T  3 

w  *n  •-»  S 

CCS 

u. 

•see 

CCS? 

er  s:  o  e 

o  c  c  c 

c  c  c 

u 

a*  lb  ul* 

ll  a.  a. 

lb  uj  lb  uj 

UJ  u:  UJ  a 

lb  ,  U-  u. 

0  «-<  c 

r?  N  <C 

NT  r*T 

T-  CM 

®J>  X 

tu 

«  <= 

TT  0  C 

N  C  T  C 

H  <0  IA 

▼  fN.  X 

X  X  r- 

N,C  CO> 

T-i  -o  -o  c 

*0  .'O  0 

X 

W  O  IT 

=  0  T* 

CSOO 

TCHJ1 

e  f* 

X 

r*  cvi  'C 

jn  cm  -o  : 

ftN-OO 

CC.  T-  to  fi 

0  >c  X 

X 

N«f  C 

<c  'O  r* 

co  fir 

«n  c  e 

v*  if>  e 

-J 

M  S'3> 

r  y  ^ 

t  r  >  t 

W  «c  CM  *0 

0  W  N. 

u. 

cm  cm  K 

CM  CMfw 

r*  ▼  •O  C 

H  O  fO  'O'' 

f.  to  CM 

z 

«  '#  i 

•  • 

•  •  «  • 

•  •  •  • 

•  •  * 

“■ 

c,  c  c 

1  ■  1 

C  <S  C 

c  e  c  c 

c  c  c  e 

ce  c 

i; 

1  J  '1 

c 

f  t 

l  It 

%*  1 

X 

*3 

< 

>■ 

c. 

* 

a 

w 

•< 

V 

- 

■J*  •  « 

•  •  •  ' 

•  •  .«  • 

•  •  • 

CCS 

s,s  c 

- 

s  s  s  s 

“  «=  =  - 

o  c*  e 

c  c  c 

c  c  c 

C  f  f  c 

C  W  f'  c. 

C  C  c 

®  c  c 

c  c  c 

cc  x  c  e 

c=,cc  e  c 

e:  c  c 

x  x  a 

LL  .  U.  L_ 

lU.au. 

ir  x  a  .a 

m  u.  u. 

in  c.  ir 

'0  a.  0 

IT  CV  tf«,CM 

Cv  C rTft 

^  C  N 

cm  rv  jv 

n  o  jr  o 

is.  ec  0  cm 

C  CM  o 

0^  0“ 

ic\  cm  <r 

cKfifi 

O  f  T  W 

CM  <5  vf\ 

*C  •”  fN 

cc  s.  a 

n.  <  it. 

f*  ▼  o«e 

c.  tt'CM 

in* 

O  Nt 

MfTIfl 

-T  10 

- 

0  0C 

O  S  « 

“ 

rl(C  OS 

t-  r  .  4^- 

X  f4  —1 

0  CV-jC* 

»  CO 

&  *  TT  <T 

C  w  Wi  - 

O  0  « 

▼  O.  W 

to  CM  *H 

■■O'w  f*  tO 

T  T^'V 

CM  rl  ri 

•  * 

*  •  • 

• 

•  •  » 

o'^;x- 

c  c  c 

e  c  o  c 

c  o  -o  c 

<C  C '  C3 

i  " 

tl’ 

1  9 

1  1 

1 

M/  + 


y 


(27) 


f  0 


f 


r 

r 

E 


[ 


[ 


j 

i 


[ 


[ 


[ 

E 


t 


o 


l 


HN  HO 

HNHO 

OC 

a  o  o  0 

0  0  0  9 

' 

Ul 

rt 

CL 

til  fit;  us  u; 

U*  UJ  U*U> 

X 

®_cm 

-< 

«c«  to  fN 

T  ^  MO 

u 

NC  CN 

*C0*>  *3  N* 

n  cm  w  ^ 

N*  *4  * 

rs.  cmp-om 

*- 

~ 

f*  «M  9-  t* 

^  o  o  «*- 

•— 

»0  CM  **  ® 

^  CM  ||N»S 

■X 

nhnh 

0',f»  N  * 

•  •  •  • 

«  •  '*  * 

>- 

o  c  c  c 

C>  C ,  C  ■  C- 

X 

1  1  1 

11  t  1 

c 

u. 

X 

— 

H  »><  H 

^  H  X 

M  N  N  a 

ec.itt  w  •= 

CM  CM  CM 

3  C  C 

C 

coco 

O'O  o  o 

O  C:  O 

Ui«L  o< 

UL.  Uy  UL 

U*  Ui  Ui  U. 

UJ  IX  lli  LL. 

lit  UL  UJ 

q: 

TMW 

«3> 

»'<c  oo  inf 

WCMN 

•— 

cm  *c  *■ 

o  *  > 

rO  *oAc  fs. 

CM  c  CM  ® 

pv  W  O 

to 

O  S\ 

=3  ▼  r4 

vOO®>CM 

3NH 

*c,  r>.  o 

ITV  C  «•< 

c  ^'On 

^-0,tc 

f-i'O  m 

•u  ~  -5  v  3 

o'm  *n 

=  Mf»  = 

®>to 

■t  o  it 

-J  C 

wac 

o.  O  CT 

NV«C 

0  !0  rH 

X 

T  ri 

»'*o  3 

lO  3  XN  3 

CM  f«f  C  O 

H  ri  H 

—  ii  *c  cm  to 

CIOT 

H  Wnlf 

HMnt 

t-i  t* 

CO 

• 

• 

• 

•  •  ••' 

•  •  •  ;• 

•  *  ,  •  "# 

•  •  • 

CO 

■< 

o 

c 

.o 

o  c  c 

•e  o  o  o 

c  e  c  o 

coo 

x 

x> 

1 

t 

•o 

III  1 

1  " 

▼ 

«  v 

•o 

1  1 

X 

u 

CL 

X. 

>- 

V>» 

>• 

>r 

>- 

X 

c 

II 

X 

tc 

X 

t£ 

w 

-4 

»■» 

n  H  n 

CM  CM  *M  O 

CM  CM  CM  3 

CM  *H  CM 

tl 

Ui 

CO 

fO 

o 

c 

o  c  o 

▼ 

3  0  0  0 

▼ 

COCO 

a  c  o 

X 

10 

X 

X 

*« 

Ul 

«— • 

II 

u: 

Ul  1U 

It 

u:  u;  (u 

II 

to  UJ  UJ  Ul 

it 

U>  Ul  Ul  UJ* 

M 

UJ  u;  m 

L5 

X 

X 

** 

*o 

r 

•O  -  .f) 

f*.  f,  CM  O' 

3  X  CM*  *0 

UN  tO  O- 

>■ 

X 

H- 

IX 

ao 

X 

a 

to  <  O' 

UJ 

*■<  cm  x  m 

Ui 

▼  ®  fi 

UJ 

▼'«  CM 

>x 

> 

X 

CA 

'M 

C 

3^ 

x 

\  Ne 

Nl 

CNXH 

NT 

*£  IN  3  3* 

v  a.c 

3 

z 

M 

;* 

CM 

3 

e  ^  o 

c.  pv  w 

MB 

OifNt  C 

•— 

O*  CM  C 

to 

-a 

U 

10 

£ 

X 

10 

w  *^  .r 

in 

£•  r^-  r-C  CM 

CO 

SHTN 

If. 

▼  CM  K 

►- 

w 

u. 

c 

X 

O' 

<r 

CVN 

CM  X  O  00 

X 

CM  ®  tO 

X 

■— 

X 

in 

-N 

“ 

■X  r-  3 

•** 

**"» 

3  r  -•  r 

»•* 

3  ▼  CM 

a 

c r 

O’ 

s 

r- 

A 

CM 

CM* 

•C  r^-T- 

n 

w  cm  ^ 

r* 

H  Pi  r! 

IT 

^  X  r-* 

MX 

< 

J- 

;• 

“• 

• 

•  *>  • 

•  •  *  • 

«  «  *  • 

•  •  • 

o 

-1 

X 

w 

S 

s 

o 

3'C  3 

w 

O  O'O  o 

»— 

O  C  3  C 

ir* 

3  0  0 

— » 

-4- 

T 

-L 

X 

f  1  1 

«L 

1  1  1 

X 

1 

u. 

II 

u 

c: 

Cl 

Cl 

Cl 

L. 

«£ 

U. 

*"* 

C 

X 

*— 

W  *-  -4 

r*  f»*  C 

T-i  -4  -*4  O 

H  H  H 

c 

w 

C 

c 

cr 

C.  C  C“ 

C-  o  o  c 

c  c.*  o.e 

c  c  c 

Ul 

\ 

U. 

u. 

u: 

Ul  Ul  U. 

U.  UL  U.  UJ 

U.  Ui  LU  Ul 

LL;  UJ  to 

*; 

r- 

vC 

•c* 

•r  O'  3 

▼  to  ^  rv 

M  ®  N  ^ 

UN  f>.  ‘S' 

z 

f*3 

■c 

> 

o  ®  o,m 

^  ®  w  r: 

MOM 

CL 

*c 

o 

< 

a:  ^  t 

o*  n  s  o 

Crt<o 

CM  tO  *-t 

X 

*” 

zz> 

'C 

o- 

o  cc  »-*  r*» 

to  to  c\yr 

*0  X-,0 

-J 

V 

.n  o>  ■? 

r-  v  o 

C 

f*  N*  r* 

u 

S\ 

>o  O  X 

W  if  ^  O' 

W  ▼  3;  If* 

O'  Mn- 

z 

S' 

nj 

CJ  «  CM 

if*  r,  xi  un 

=  x  r  f> 

rO  X  CM 

r. 

If 

0“ 

^  rv 

U*  O'  Mr* 

If  X  Mri 

r>  M  n 

• 

• 

• 

•  •  <•  • 

•  •  i  • 

•  •  6 

— 

1 

o 

coo 

1  1 

COCO 

1 

COCO 

1 

COO 

1  1 

0 


fl 


I 


I 


s 


*>  < 


V. 


* 


U  C'-(.  u  l; 


0 ) 

,  /  .  .  C. 

2.5  PBOStAM  LISTING 


*  '  FORTRAN  DECSf  1 

CHAIN  PROGRAH  STRIP  -  AERODYNAMIC  INFLUENCE  COEFFICIENTS  BY  2 

INCOHPRFSSIBLF  STRIP  THEORY.  *  3 

HITH  OR  WITHOUT  CAMBER  4 

WITH  OR  WITHOUT  A  CONTROL  SURFACE  „  5 

NCAH  =  0>  CAMBER  NOT  CONSIDERED  NCAM  =  1,  CAHBE li  INCLUDED  6 

71c,.  7.2,  Z3  ARE  PERCENT  CHORDS  OF  C.F.  1*2*3  RESP.  ON  EACH  STRIP.  7 

r  A  NEGATIVE  1.7  INDICATES  presence  of  a  control  surface  on  the  STRIP  8 

C  9 

n  I  Mt  ns  I  ON  TITlF(24),Xl(25),X2(25),X3(25),Zl(25),Z2i25;>,Z3C25),  1  u 

.1  a  m  ( 3 ,  .3 ,25  ) ,  »M  ( 3, 6  JrHU(  ?5  ) ,  B  ( 25  >  ,  SCALER  ( 25  >  *  CHi<4»8*  25  )  *  PARTA  (4,8),  13 

X  lA(25),C(.?5).f  J0X(3fi),AHT(3,3*25),Df:L7Ar{25),ri(25X,EY0X(7),  12 

.1  f URl  ( 25 ) ,  Uri  (25  >  ,-TA  12  (25 ) ,  TBP3  ( 25  )  *  Tfc  12  (  25 ) ,  THR2 (25  ) ,  TAR2  ( 25 )  *  3  3 

4  <  HPl  ( 25  )',  TRIM  ( 75  )  ,  TO  1 3  ( 25  )  ,  T H 1 1  (  251 ,  T A  1 1  (25  )  i  TBR2  f'25  )  ,  TB 1 1<  25 ) ,  14 

5-1  R«1  (25), El  0«3(25)  .T:LBI2(25 ),EHBR1  (25 )., EMRI3 (.25 )  ,ELRR2(25),  18 

*  6-T  n  1 1  ( 25  ) ,  Fi  r  i  3  ( 25  )  >  FMRR2  ( 25  )  ,  EKR  ( 511 )  ,  Z2*( 25  j *  XLt;25  ) ,  XH  ( 25  )  ,  XT ( 25 )  16 

7  .  HI!  ( 25  ) ,  HL2  ( 25  ) ,  Hl.3  (25  )  ,  MT1  (  25  )  ,  HT2  ( 25  )  ,  HT3425  >,CAH  (4,4, 25  ) ,  1 7 

H‘10H(  25 ) ,  Hill  (P1' ) ,  HH2(?5  ).,  HH3 (25  )  * CA,MT  (4 ,4 , 25  )  ,  XNl  f  25 )  >  XN2  ( 25 ) »,  1H 

.  9*N3(26),XN4.(25),XM5(?5),T1(25),  72 (25  )  , T3 (25 3 »  74  (,?5 )  i  T5 ( 25 3  ,  CM(4, 8 )  19 

f  -  2.5 

1  f ORHAr(4l4/5F12.8)  24 

2  FORMA  T ( 1H11  25 

3  f*ORHA T ( 5F12 .8  )  30 

4  FORMA T(lH0  3 9 X , 81 H  AERODYNAMIC  INFLUENCE  CUEFFICIENIS  BY  INCOMPRES  35 

1SIRI.F  STRIP  THFORY  HJIHMIT  CAMBER/// )  36 

5  FORMAT (1  HO  54 X#  1.1  H;  INPUT  DATA  //1H  4 8 X ,  12,7/4  STRIpS/lH  48X.12,  40 

1?2H  RFIlUCEn  VELOCITIES  AS  HO  45X, 1 5HC0S 1 NF  I  AMBDA  1E14.6/1H  38  41 

2S22HRFFERFNCF  SEMlrCUORD  =  1E14.6/1H  49X,llHSE.Ml-SPAN  =  1E14.6/  42 

3 1 U  46X,  14HSUKF  ACE  AREA  =  1E14.6/1H  53X,7HC  /jA'R  =  1E14.6//3  43 

6  FORMAT ( 1  4Hfl  1 /K < R  )  =  2E1 8 .5, 4E1 5.,5/T  14X,2E18.5, 4E15 .5  )  >  45 

7  rOPMA  T  ( 6E12  i:8  3.  46 

ft  KIRHAUTMO  4X,9HSTRIP  NU. > 7X» 11HDEI  IA  Y  (  I  ) , * lXi 4HB ( I  ) ,  50 

1  l  0X,5H?1  ( I  ),1!IX,5HZ?(  i),10X,5HZ3(  I ) ,  1  0X» 4HD (I  ),10X,5hCA(I  )//(  19,  51 

2-23. 5,  FllJ.5, 51  15.5)3  52 

„  9  ‘ORMATMHfl  48X,8H1/N(R)  =  E13.6)  55 

m  *  0  P  M  A  T  ( 3  H  0  53X. UNSTEADY  C4SE//43X,  33hl/K  ( R;)  =  LARGE  NUMBER  (OMEGA  ,60 

1  =  h.))  '  61 

11  “IJPMA  T  (i  H  29  V  ,  4F1  8.8  )  65 

12  FORMAT (12 A6)  66 

!  3  FORMA  I  M  HO  48Y,  19HNIIM8EH  OF  STRIPS  =  I  2  »  *  7U 

1  4<  F nRf"A  f  (  12h  2  E 1 6 . 8 , 1’X ,  2E10.8,  I  X,  2E16.DT  75 

1 5  rORPA  f  ( 53110  .  CH(J?o  80 

1PH)  SIZE  =  II, 4H  BY  II//)  PI 

16  rOPMAT ( 1  HI  12A6//1X,12A6///)  86 

85  lOPMATTlH  ?F  I  6 . 8  ,T  X ,  2F 16 . 8 ,  l X,  2F16 . 8 , 1X>  2E16 . 8)  87 

216  formATC  1  HI)  7  X  " ,  7 8 M  AERODYNAMIC  INFLUENCE  COEFFICIENTS  BY  INCOMPRES  91 

1SIRIF  STRIP  THFORY  Wl fH  CAMBER/// )  92 

220  iORmAICIHO  4X,9HSTRIP  M0.,7X,11H1)EI  TA  Y/  (I  ) » 1  LX,  4HB  (  I  ) ,  10  X,  5HZ1  ( I  )  98 

1  . 1  n  x,  5H72  (  I) .  TO  X,  5H73LI  ),17IX,5HCAH  ) // (  1 9,  F23 .5,  E18 . 5, 4E15 . 5  ) )  99 

r  •■T-*'.  INPUT  1X6  I A  AND  PRINT  100 

;>n  HP.  AH:  (6. 12  )  (I  i  IT F,(3  h,  1=1  i,  24).  101 

* P  I  IF(6,1A)(  rm.F(T'Y,T=i,241;  J  0 •> 

,>l: A P  ( 8 . 1  )  NF  A-:.  I SZ,  JS7 ,  NOPlih J»  COSLMI) ,  UR,  S»  CAPS,  CBAR  1  ? U 

HF-A|5(5.3)(nu  I  A  Y,(T  >,B(I  ),,ZKI  )„/?(!  ),  iZ(  I  ),  l  =  l,„ISZ)  335 

•>FA«i(5.7)  (EMM  D,  1=1  ,  JRZ)  140 

•10  50  1=1,1.5  7  143 

50  ‘>2A(  L)  =  ABS(  17  <  I  ) )  142 

IF(MCAM.PU.Io)  GO  10  300  145 

■JR  I  I  F  { 6 , 4  )  15  0 


(29) 


00  TO  310 
3Un  MRHF(^216) 

310  RR 1 TF ( 6, 5 )  T  SZ » JSZ* CQSLMD, BR,  S»,CAPS»CBAR 
CAKRFR  OR  Nd  CAMBER  OPTIONS; 

IFCNCAM.EO.l V  00  TO  500 
00  55  1=1*1 S7 

0(1  )=2.0«B(!n*(Z2ACF)-ZTC  1>> 

I F  ( Z?  (.1  )  •LT«0 . 0  )  GO  TO  37 

cam)=o.o 

00  TO  55 

17  0 A ( I )=?.fl*B( I )*(1.0-Z2A( 1) ) 

55  C0N1INUE 
00  TO  550 
5tm  »io  58  1=1, ]S7 

<1  (  I  )=  2.n#R(T  )*Z1(  I?) 

V? (|)=  2.0*BM  )*Z2ACtf  ) 

<3(-n=  2.o*B(>n*z.ic  r<V 
xi  ( i  )=n.n 
y}i{i)=  r(I) 

XT( 1 )=  2 • 0*B (  I  3 
I  F  ( 7.2  (  I).LT.n.n  )  GO  to  51 
CAM  >  =  ft.,0 
0(1  TO  58 

51  TAM  >=2.0*B(  f)*<1 -0-Z3(  I  )  > 

4H  n  F  ( 6 » ?  2  0  >  ( .  I  .  DELTA  Y(  U)  ,  B(  I  )  ,  Z1  <  I  > »  Z2A  (  I.) »  Z3(  I  )»CAC  1  )t  1=3 

55  0  <>R  ITFC6, 0)  (  l  .hFLTaY  ( I  ) ,  8  ( I ) »  Z1  ( I) »  2?  A  <i  ) »  Z3  (I  ) ,  D<  1 ) ,  CA  (I  )j 
551,  mitF(5»6)  (FtR(J).  J=1*JS;Z) 
in  ca=o 

Hi  ?\  1=1.157 

i'F  (  A  F AM . Ell «  3 -)  00  TO  510 
iF  =(  C A t  I  )  )  21,22.21 

510  I  F  CO  A  (  1  >  >51.1  .51  2,5li 
2V  iti(I)=3 

•‘0  CA  =  1 
NOTH  23 

511  'l II (  |  )  s  4 
>'ior:A=l 

o  n  LO  23 
51?  -Jll f  1  >  =3 
on  TO  23 

22  N 1 1  ( 1  )  =  2 

23  CONI l NOP 

»I  =  3.141592-55 
nill  =  1.0/ FI 
l  j  I  ?  -  P(1  l  »*2 
IF(NCAM.fii).1  )  00  TO  finn 

c;ca  con  =(pu  cosimri/u  hr**2>  )  *s> 

•in  26  1=1,157 

Of'.  A I  FR  (  I  )  =  5C.A  CON'  »  (8(1)  *«2  >*  DILTA  Y(I) 

m  p  =  !i ( !  >/n<  i ) 

AH(1  ,  1  ,  1  )=  1  .  ’I 

\M(1  ,2,  I  )  =-rl  il'«Hll) 

am  (2,i*l)  =  n . n 

am  (2,2.1)  =m n 

IF  (  CAW)  )  24,25,24 

24  AM  (d,  3,  I)  =01  n 
WCft  =  H  (  I  )  /  C  A  (  I  ) 

.aM  (2,o,I)  =-Hl  D-H1F.A5 
AM  (3,1,1)  =1.0 


,  ISZ) 
1=1, ISZ) 


( 


<3\ 


AM  (3,2,11  =  0.0  425 

AM  (3,3,11  =aiTA  _  *  „-430 

25  Ms  ,NU ( 1  )  43 5 

HO  ?.f>  K=1,N  44  ii 

on  26  ‘1=1,  N  445 

26  AMT  (K,T.»I1=  aM(L,X,  I  1*SCALER(U  450; 

GO  TO  41  459 

600  nAMnnN=(Pl*.(:h«?LKD)/c(  <BR«*2)*S)  460 

in  610  1  =  1,1  *?2  465 

GONl  I.)=nAMGOM*(H(  1 1.**?l*nE:LlAY(  I)  470 

-u.1  ( i  )=(xlc  i  )irx2rm*(xi  < D-X3(  i  >  i/uxki )-x?t  m*(xid i-x3u;i n  -*75 

HT1  (  I  )=  (  XT ( 1  1-X2T 1-') .) *  ( XI  (1  1-X3(  I  1 )/(  (XI  (  I  )-X2  Ml  )'*'(  Xl (  I  >-X3(I)ll  4,40 

BL2  M  1=  (  XI  (  I  }— XI  ( I  11  *  (XL  ( t  1 - X3  ( l/)  )  /  ( (7  2  (  I  )  ^  XI  ( XI )  *  (5X2  (  I  X- X3  ( I  11  )  A  B5 

MT,?(  I  )={XTX  n-Xl  (  II  j * ( XT  (  I  1-X3X  I  )  )/(  (X24*M-X1  M  )  )*(X2(  I7~X3(  I  ) ) )  4911 

41:3  (  11=  (  XL  (  1  )-Xl(  ]  1  ) *  ( XL  (  n -  X2  ( I  ) ) /  ( ( X3  ( ,f  >-Xl  ( 14  ) *1  X3  ( I  1-X2  ( Ill  )  4 95 

‘IT3  ( I  )  =  ( XT  n  )  -  Xl  ( 1 ) ) » (  XT  ( 1  )  -  X2  (3  14  /  ( ( X  3  ( I) -X1M7 1 *  C  X  3 (1  i  - X2 (I > 1 )  5  0  u 

•IH1X  U=(XH(  1  )-X2(  I  )  )«(XH(  I)-X3(l)=)/(  (Xl(  T  )-X?;(  n^*(Xl(  I  )-X3(  I  11)  505 

4H?(1  >=(XH(  H-Xl  (  I  ) )  *  ( X<H  U'l-X3  (T1 )  /.  ( ( X2.  ( I4-Xi  {T1;}*XX2T  I  )-X3(  1)  )')  510 

4H3  (  I  )=(XH(  I  1-Xld  )  )*(:XH(!  >-X2(  ])  >/(  (-X3(  I  )-Xl  H  I  1*1  X3-q  >,~X2(  IT)  )  515 

CJHI4  /1 ,1  )=i1A»Hl.1  (>1,)+.7S*HU(X<j  '  520< 

G AM ( 1 , 2 , I  )=.75*H|  2(  1  )+.2t«HT5(1i  )  525 

GAMil  ,3,1  )=.75*HL3(  ri  +  .25*HT3(  II  530 

G  AM  ( 2 ,1  ,  U=.5*(HT1(  I  .),-HL1?( !:)  )  535 

•:am(2,7,  U=.5*(HT2( i  )-Hl.2(i  ) )  540 

GAM (2, 3,  1  i=.5#(HT3(  !4-HL3< I  )  )  546 

GAfU^-l-disHHl  (  1 1-.6f  (HLKiH-HTK  I.)  )  650 

GAM ('■3  ,^ ,  I  )=L:H?  (  H--.5-* (‘HL'2  (1  1+HT2 ( I ) )  556, 

GAM(3,3,  1  )=LH3(  r)-.T»*(HL3(r.H-Ht3U:i  )•  660 

>  F  ( G A  ( I  )  16 0  ?/,  6  05,-60?  -  565 

602  GAM (1,4 U  1=0.  IJ  570 

GAM (2, 4,  I  }=0„0  575 

GAM (3, '4,  I  )=0 . 11  5 flu 

GA‘M(4,Xi  I  )=-GU  )*(X3(4.)-X2(  I  1  )/(  (XI  (1)-X2<  1  ).)«(X1  (  I  )—  X3'( .1  11.1  586 

GAM  (4, 2,  1  )  =  -RU  )*(,Xo(  I  l-Xl(l’)  )/( CX?(  I  J-.X1  (I  ),!#,(  X2  ( I  1-X3  (  I  1 1.)  590 

GAM (4, 3,1  1--H (  I  1* (2 .,0  #X3 (  1  1-Xli  I1-X2 (  L> )/( (X3(I 1-XlU  1)*(X3(  I  )-X2(  595 

1  n  )  )s>?(  1  )./ca<  1 )  60  0' 

g  a M  ( 4 , 4 ,  T  )= n  ( .1  >./ n  a  (  I )  hn5 

6115  M=Mil(M  fi^ll 

'HI  610  K=T.»  N  615 

•’in  610  l  =1  ,  fl  620 

6in  GAM1  (K,l, »  1  )  =  GAM(L,K,:I  )*r,ON<  ]  )  625 

4.1  IF  (NOG  A  )  31  ,3.1,27  63  0 

27  <F  (  ,JS2-1  )  29,20,29  635 

2B  1 r  (  FKR(l)  )  30,29,29  640 

29-  GAI  I  ORANGF  645 

i  (  I  S7 ,  G A »  M ,  1  w»<i  ,  THR2, 1 H 1 1, TA R1 ,  T  AR2  ,  TA U  ,TA  1 2,  1  BRl  ,  TWR2 ,7RR3 ,  650- 

2IRR4,  Toll,  T'I2,  T0I3,  >FLRH1  >  FLBR2,  fcLHR3,  FLBIl,  fcLB 1 2 ,  ELS, 1 3 ,  655 

3FMRR1,  RMHR2,  FMRil  .),  660 

3  0  I  F  (  U  G.A  M  )  3 1 , 3 1  ,.51 5  6  6,1 

515  GAM  GST AM  ( 1 47 ,  CA,  P»  XN1 ,  XN2 , XN3 , XM4 , X N5 , d  »T2 »T3, T4 ,35  1  662 

31  rONITNOF  666 

iSFO=0  666 

on  70  J=l,  IS7  670 

MO  60  1=1,  1517  675 

1-1  K=FK'R(  J)*WR/R(  1  )  68  0 

IF  ( F  K  I?  ( J )  )'  37,36,34  685 

37  !F(IMGAM.FlJ,1  )  00  TO  9f  686 

GO  10  90  687 

33  1=0.5  690 

0=0.0  69,5 
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GOTO  35 

34  PK  now=i  ,  o/m  K 

"All  •  RFSSEI  ( FK  NOW,  1,  EJflX,  FYflX,  0) 

'  =  (FJnx(2)  *(F.JOX(?)  ♦JYOX(l))  +  EYOX  (2 )*  (E.YOX  (2)  -EJOX(l))) 
1/  (.(  F.IOX  { 2  >  «-»*  YHX  t  i  >  >  ##2.0+(  EY0X<2>  -EJOX  (I) >**2,0) 
>=-1.0»(FYOX(?)*EYOX(l)+EJOX{2)#  EJ0X(1) )  /<  ( E JOX (2 ) +EYUX ( 1 ) 

1  i  »*?.»  -»  ( F  Y  OX  (  2  )  -EJOX  (D)  ##2.0  ) 

35  r-ONTINUE  ' 

UK?=F1K*E1K 

•< M (-.1  ,  J  )=1 .0  +  ?.0*R#E1K 
•»H(.J  ,2)s-2.U*F*fclK 
*»M(1  ,3)  =  0.5*?.C,*fi*E1  K-2  .,0  *F*EK2 
HM(1  ,  4  )=-FlKr?.0#F#,FlK-2V0#G#EK2 
=  I)  .  «=. 

=n.u 

•IM(2,3)  =0.1/5 

»M(2„  4  )  =  -£U; 

>nrA(f))  3/-.  5411,36- 

36  -'MI  J  ,5  )=l'.*EJ=ff#FL;RI<l  (  I ) -F»FK2#EU8R2 (  I  KELRit.K  J.) 

»M(i,6)=-F*llK*fcl.  Gil  (I  )-G#FK2*FLRI2(a  )  >EtK»FL  8 1 3  ( I ) 

•«M(7,5)  =  EHBh1  <  I  )-EK?*FHBR2(!  ) 

<M(*2,6  )=-FlK#l*Mrtir(  I  ) 

l )  =  T ilHl  (  I  1  ♦R#FJ,K»T»IR2C  I  ) 

I  #1  1»  *  rHjl'f  I  )' 

'Hf 3,3)  =  TAPI  (  I  )*!(G»I  1  K-E#EK2)#TAR2<  1  ) 
i M  ( 3 ,4  ) =  (  -  F  »  P 1  K  - R  # F  K  2  )  «T A  1 1  (  I )  - F 1  K«-J  A  J  ?  (  1  ) 
iM 1 ,1  »-5  )  =R#F j  K e TiIRl  (  I  )-F#FK2#TMR2(T  )^T8R3(  I  ) 

1  -FK?»TBk4.(M  ) 

■  Mf  .1,6  )  =  rt  *i  1  ••  *  foil  (  !  )-:0*FK2*TBI.2.(<? ) -I- IK#  TB 13  <  I  ) 

54(1  •(  (^rAM.FO.  1  )  ft'»  If!  541 
ill  Tii  40 

541  -Mf 1  ,  P):PM(|  ,1  ) 

•  H(  i  ,  [>  )=r»n(  i  .  ?  > 

•H(  1  .  i)  =8*1(1  ) 

»-M(  I  .  4  )=Brt(  J  .4  1 
f,M  (  ?  ,1  )  =8  M  ( ?• .  I  ) 

’M( V,2 >=8M(?,2' 

:mI7, 3  )=8M(2.  l) 

•M02.4  )=HH(2.i  ) 

"M41 , 5  )  =  .  75  ♦<*!  K*rji  2.  (1«FK2»F 
• - M'(vi>.  6  )  =  - 1  iNl  ♦  2 .  ii # h  K 2 * 0 
GMf.Rx,  5)  =  .J7‘»h>K’X 

•  M  (  2 , 6  )  = .  5  #  1-1  K 

)  =  ./5  +  Fir«0 
"M  r .,5 )  S-M  f.  «t 
•M(3,3  )  =  .37‘/  +  i  1  K  *i;-t-Kv'*l 
.•Ml  3 , 4  )  =  -Flr  -  I  1  K#f  -V  K?#i> 

M  ( 3 , 5  )  = . 58 J  l3333> . 5»E1K*R+.5#FK?*FK?*F 
:M(3,  6  )  =  - . 5 3 F  I  K *F  ►  F  K  2  »  R 
I >F(I  A  (  I  )  >84>’ .  40 , 542 
54?  >:M(  1  ,  7  )  =RM <  1,‘>) 
i’M(  l  ,i\  )=H«(  J  ,6  ) 

•:m ( ?,  7  )  =8H( ? .  5  ) 

GM(?,  8)=BB>(?  .  6  ) 

•M<3, 7  )=“F  #1  K2»XM1  (  J  >+R»ElK#XN?(Y  )-f-K?*XN3<  I  )  +  XN5<  I  ) 

’M  ( 3  >  8  )  =  -G*l,KV»>XN  I  (  |  )  -  F  »F  i.K*XN2  (  I  )-tlh*XM4(  J, ) 

’MM,  I  )  =83(3,  I  ) 

•  ’MM-,  ?  )=8i-l(3,2) 

■M(4,3)=HM(3,  U 
i’m  (  4 , 4  )  =«M(  J  .  4  > 

.’MM,5)=-I  *1  K5«  fli  |5j):+i;#Fl'K#12'(  T')-EK2#T3  (  I  )  +  15  (  f  ) 
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4,6  )=-i;*r-KP*Tl  ( I  )-F*FiK*72  ( T  )  99b 

CM{ 4 » 7 )  =Bf1  (3  ,5  )<  109° 

ON ( 4 , 8  3  =BH ( 3,6) 

mi  in  4n  '  lntl6 

STFAOY  CASE  OPTION  5nn9 

pn  .'1  =  NII(  I  ) 

=  1U2U' 

no  38  JJ=3,N?  '  1H25 

•ill  38  u=i,?(  J-f'3" 

38  ^M(I I, jj)=o.n  3  03b 

91  '0P|:-{2.IJ?S«C»AR/{:AI’S)*BR*BH/(8(  JT«B(  I  )  >  103ft 

:iM(1,l)=T2.ll*CfrRR  104b 

4  F  (A' A  (  1  ) )  39  >  9?,  39  1  (>5U 

39  •«:'(  I)=1  .IItCAU  )/i1(  I  )  10611 

,JH-ARCSC-C.(  I  )  >  106b 

rn8P=-r.(  u  1070 

;inp=sjm(ph)  .  .  io7b 

»M(3,3  )=-ri  J 1  «  ( (  P  i-l,H)*(-l.fl  +  COSP  +  nOSF  l+4  3.n-n0SP)«STNP)«C0RR  1118b 

~<H(1  ,b)=-H.lt>i'l  M«(PI-PH  +  SINP)'«r.OKP  "  1090 

■»m(p,5)=-piji«sinp»(1  .n-nnsp)*C0RK  in9b 

«j)=RH<3.3S«Pll1f  (PI-PH*SINP)-rP.I3  HOU 

1  -  *>  I  N  P  *  ( 1  .0-.i:nSP)#(PI-PH-SINP)*n(lRH  310b 

9?  iMNCAH.EU.flufl  T1I  93  13  Oft 

•in  -in  4ii  no/ 

93  J  =  HII(I)  110M 

i?  =  b*?  13  09 

.HI  94  JJ=1,-h»  1110 

•  ()  94  11=1  ,f-  HU 

94  'Mill  ,  J-* )  =0-.  0  -  3.11 2. 

"fKI  ,3)=UM(1  .  '  >  1 113 

#MU  ,5)=2i  0#ni'RR  1  33  4, 

) =cuhk  m& 

fMf.V,  3  )  =  -!..  ji*  riOKH  1116 

•■K(3,  b  )=CORr  *1  .*>  '  113  7 

•  F  ( I'A  (  I  )  )’9b ,  4  :• ,  9b  1118 

98  "Mf  3  #  7)=PMf1  1119 

**W2#7)=RM(2,v)  ‘  1 T?  0 

ff-l  (4  > .V)  =  R|1  (3.3)  3  3  21 

•  ‘ 7)  =011(3.'.)  11  ?2 

•:m(43,7  )=COR|.  s  (  -,XN1  -  SN3.)  1123 

(  4 , b  )  =C OR1-  «• .  -  1,1  -  f  3  )  1124 

4  0  *:flMi'!NIII  •  13.2b 

I  F  n-r.AM.Fl).  1  )hfl  TO  r- flu  112h 

;fmi pa.ti-  ah;  *a trici-s,.  print.  results  and  punch  output.  n?/ 

■•All'  WMAim  (  A-'.OM,,AI'T,'Cni  ,  NU,  1,3,  ft)  1329 

10  10  All  13  30 
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if  o  m  h  in  ri38 

80  0  HR!  T  F  (  ft  ><2  lift  >  -  1139 

Hl!1,  M'U'KR(J))  61,62, ft?  '  1140 

•61  RUF(6,1il)  114b 

TO  10  63  1150 

6?  *,’R  I  T F  (  6 , 9  )  l-K  ^  (  I )  1155 

63  wW  I  1VE\(  6 , 1  30  1^7  116U 

in  68  >1=1,19/  116b 

'BIlFdl.fi)  I  .  MH(  I  ) ,  NiJ  C  1  )  1170 
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■  1 O  r-fj  *  ?  '  11  *1 II 

IFU-KR(J))  64,66,66  1  1.H*# 

64  mH  65  K  =  1 , N  1  ! V n 

65  JRlIFlfi.ll  }  ('MUtK.l  .  [  >,L=1,N?,?)  1 175 
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66  if)  67  K  =  L,I  |?l!5, 

«ff  NCAf'.tU.  1)  011  fn  t8n  K'lJI 
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60  iscii=nSI-0»1  1  ?4 1 

•All  iMINJ  (I  K  '(  J)  *K(  »CM1,  iS7-,NSrO)  I  <’  1  ■* 

7 n  '-nMimif  '  .1  ?5 ii 
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B*V  )*(n.a>*4.(l*  CUS2)+  (Pl-PH) 
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2ft.il 

2  ll  1  "r 

2in: 
?i  i ; 
2 1. 2 1 
Pnpi, 
2 1  .<  ii 
2n,la, 
2H4it 
2  ft  Ah 
2  ll  S  n 

2 1<  >i  1/ 

2  ii  A  ii 
PllA*^ 
?07ii 

?n?f- 

?(lRji 
2('Htj 
2  0  V  H 
?H/b 
2i  mi 

21  (1*5 

?11fi 
211  ij 
21?n 
21=2‘; 
2  I  3  n 
2]  3,5 
2  1  4  l) 
2 1.4  '•> 

21*)  ii 

2.1  M, 

2.1  611 
2165 
?1>7  ll 

21  /  ‘i 

21  HJi 
21 H  M 
2  1  9jl 
2 !  95 
?2ll,(i 
22H  S 
221,0 
2,2 1  r> 
22211 
???*> 
2  23  (>’ 
2  23a 
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S  FORTRAN  DECK  3001 

CHMATHV  COMPUTES  AERODYNAMIC  MATRIX  BY  PARTITIONS  300? 

SUBROUTINE  W  HATH  1 ( AH,  BN# AMT# CHI # NU#  l#Hl»M2)  3005 

DIMENSION  AH < HI. HI# 25 )» BM< Ml# M2), CHI (4. 6. 25), AMT (HI# Ml* 25)  3010 

1 ,PARTA(4»8),NUJ l)  3011 

NU1  =NIJ(T>  3015 

MU?  =  NtV( I)  *2  3020 

HO  80  K=i#NIM  3025  ' 

HO  80  L=1 # NU?  3030 

PAPT  A(K,L)  =0.0  3035 

HO  8C  J=l#N1i1  3040 

80  PART  A(K,  L)  =  PART  A ( K ,  L ) <•  AH ( K,  J .  1  )  •  BHCJ.L)  3045' 

no  83  K=1,NU1  ■'3050 

■Hi  81  L°=l, NU? , 2  3055 

4=  C  £  *3  >/2.  3080 

CHI  (K»l.»  I  )  =  n.n  3065 

no  81  J  =  1 # N01  *  3070 

8)  CHI  (X,l  ,  I  )  =CHl(K,l  .  I,)>  +PART  A  (K.2«J-1)  *AMT  (J,M,I)  3075; 

no  83=  l.=2»NU?.2  30RU 

m=  Xf>  3085 

CHI  U.l.l)  =0.0  3090 

•10  83  Jsl.NVI  3095 

83  r.Hl  <K,l»i>  =CHl  (K,L,I  )  +PART  A  (K,.2«J)  ♦  AMT  <J#H#i),  3100 

PETUPN  3305 

END  3110 


%  FORTRAN  DFCK 

f  AKCC  rr.MPUTFS  THE  ARCfOSlNF  ,  IN  RADIANS,  OF  A  pFAL  ARGUMENT, 

r 

FOVrTIOM  ARCS ( X  } 

F(Y)  =  1. 57n7g633  -  ,21460384  #  X  +  .0890456';  *  X#«2  -  .05072733 

1  *  X«*3  t  .03313246  •  V*#4  -  .021,99838  *  X«*5  +  .01261235  *  X*»6 

2  -  .rOA9l,7P6  *  X # # 7  ♦  .00096128  •  X*#8 
AX  =  ARS<X) 

IF  (AX-1  .0)  *»d  *  *50 ,4n 
<1  •*  WRITF.  (6.45=)  X 

45  FORMAT  ( /'//7H«»**»f*/26H*ERF0Rf  IN  SIJBR.  ARCS,  X  =  ,E11.4, 

1  3  3  H  liHOSF  ARS  VM.UF  IS  LARGFR  T‘~AN  177H« *##«»#// / ) 

ARCS  =  0. 
no  in  60 

50  IF  (kY)  87,, 55. 55 

67  Aprs  =  3 . 141 59266  -  SORT (1. ft -AX ) «F ( AX ) 

RO  T P  '60 

65  ARfS  -  SORT ( 1 . 0-X ) *F  FX  V 
60  RFMiPN 
F  Nil 


ARCSOOOO  * 

ARCsnnoi 

ARCS0002  5 

ARCS0003 

ARCS0004 

APCS0005 

APCS0006  < 

ARCS00C7  ' 

APCS0008  % 

ARCS0  009 

ARCS00T0 

ARCS0021 

ARCS0012 

ARCS0d.1,3 

rARCSO  01 4 

ARCS0015 

ARCS0016  ’> 

ARCS0017 

ARCS0018 

ARCS0019 
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J  FORTRAN  DECK 

CRESSF!  COMPUTES  BESStL  FUNCTIONS  (1)  OF  THE  FIRST  KIND  (  JN(X>  )  6002 

C  AND/OR  (2)  OF  THE  SECOND  KIND  (  YNU)  >  6003 

C  6006 

C  y  =  ARGUMENT  -N  *  ORDER  <0,1»2,3,4,  OR  5  )  6007 

C  FJ  =  J  ANSWERS  T  =  41  ,  COMPUTE  ONLY  Y  6008 

C  rY  =  Y  ANSWERS  =  0  ,  COMPUTE  DOTH  Y  AND  J  '  6009 

r  5-1  -,  COMPUTE  ONLY  J  6010 

c  j  6  mi 

r.  '  USFS  FUNCTIONS  A=BJ0(X)  OR  B  =  BYO;<X)  fOR  ORDER  0  6012. 

C  <AND  A=B.I1(X)  OR  B=BYKX>  FOR  ORDER  1  6013 

C  601 

SUBROUTINE  OF-SSFL  (  X,  N,  FJ,  FYV  T  )  602(1 

OIKFNSIpN  Fi(l),  FYU)  6025 

C  '  603d 

r  ALWAYS  FIND  7ERII  ORDER  .VALUES  603$ 

**  .  .  604  (i> 

l'J(l-)=BJ0(X)  604-5 

FY  ( 1)  =B  YO  (  X  ):  6  046 

1 F  (  N  )  5  ft  ,50.*  Ill  •  6  0  5it 

10  «  J(P)sRJl(X)  6055 

"Y ( 2 ) =B  Yl ( X  )  6056 

■F  (  N-l  )  50.50.12  6060 

1?  !T  .(  T  ')  14,1 4, *4  6Cfj$ 

34  >Y(3>)  =  ,2.>*LY(2)/X  -  FYT-1  )  6070 

16  IF  (  T  )  17,17,18  -  -  6075 

17  rjf.d  s  2v*l  J(2)/X  -  F!;l(l )  6080 

18  IF  (  N-2  )  50/50,20'  -  6085 

20  !F  (  T  )  24.22,2?  6090 

??  F Y (4i)  =  (8./(X«X>-1.)*FY(2)  -  4.*TY(1  )AX-  6095 

24  IT  (  1  )  24,26,28  610.0 

26  rj(4:)  ( 8 .  /  (  X*X  )-l  .  )  «FJ  ( 2 )  -  4.*FJM)/A  6105 

28  IF  (  N-3  )  50 ,  50,,.30  6110 

30  v  s  d.  -  24.-/(X*X) )  •  6115 

7  ■=  8.*(6./(XnX)-l.  r/t  6120 

I  F  (  7  )  ,  .5? .  32  6125' 

3?  fy<5)  =  Y*FY  d  )  ♦  Z*FY 62 )  '  613li 

34  IF  {  T  )  36  6,38  6235 

36  F J (5  )  =  Y#BJd)  4  Z*F J ( 2  )  614(1 

38  IF  (  N-4  )  50,50,40  6145 

40  Y  =  12.4(1.-16. /  LX *  X  ) .)  /  X  63  50 

7  =  (1-.r72./(X*XL+3S4./(X##4))  6155 

IF  (  7  )  44.4?:,,42;  616  0 

42  !- Y  (6  )  =  Y*F)U)  4  Z#FY<2)  6165 

4  4  IF  (  T  )  46,  .4  6,50  617  U 

46  4.1(6)  =  Y*FJ(1.°)  +  7«F J ( 2 )  6175 

50  4L7URN  6180 

END  6185 


«  n 


-  s  .  (/d) 


f  FORTRAN  DECK 

CfSCAH  DETERMINES  THE  CONSTANTS  INVARIANT  NltH  REDUCED  VELOCITY  5802 

INVOLVE^  IN  THE  OSCILLATORY  COEFFICIENTS  FOR  THE  CONTROL  5883 

SURFACE  UNDERGOING  CHANGES  IN  CAMBER.  5084 

SUBROUTINE  CSCAM  (  I S Z > C A *.R # XN1  > XN2 p  X  N3  *  XN4 *XN5» T1  > T2 , T3> T4 # T5‘)  5018 

DIMENSION  B ( 1 j .CA(1 ) »  XN1 (1 ) , XN2 ( 1 )>  XN3 (1 ) $ XN4 (1) * XN5( 1 ) #T1(1)  5028 

1  T3(l),  T4(l  ).T5(1  ).C(25)  50?i= 

PI=3. 14159265  5025 

R]n=l.«/PI  5038 

on  yoo  i=l,  isy  5035 

1 F  t  CA  ( I  1)706 ,700*706  5040 

706  0(  I  }  =  1 .0-CA<  I  )/BI  !  )  .5045 

PH  =  ARCS(-C( i  )  )  -  5050 

PH?=PH/2.«  5055 

•>H3=PH/3.H  5060 

>»H3hspH/.30.n  .  5065 

XN1  ( I  )  =  P1,11*(P1-PH*S1N(PH) ) .  *  5070 

XR?(  I )  =P  I&1>»  (IIP I  -PH )  /?  .  0  >  *  <  1 .  O  *2 . 0«CGS  ( PH)  )+S!N(PH2)*C2.0*COS(PH)  5075 

1  > )  ‘  5080 

VN3M  )  =  (?.0/C  i.ll*Pl>)*((SIN(PH))««3>  5085 

<*•4  t  t  ):  =PI11*<PI-PH*1.333*SIN(PH>*C0S(PH)-(SIN(PH)*({C0S(  PH3  )  )*»3 )  5090 

1  ' )  5095 

<N5( I  )=PIU*( .75*(PI-PH)#COS(PH)+SIN(PH2)*(1.0*((COS(PH2) )**2>)“<{  5100 

1 S I N ( PH30 )  )**5  ) )  5105 

M<(  I  ).=  PI  il,»  ( (  Rf-PH )«  M  •  0~2.0.#COS(PH)  )  +  S  lN(PH)*(COS(  PH)-2. 0  >  )  5110 

IP(  I  )=-U<  I  J/p.n  5115 

13 ( I )=PIll*(-rQS{PH)*siN(PH>-(PlrPH)-M  .6667*( (SINIPH) )*#3>)>  5120 

r  4  <  1  )  =  P 1 1 1 «  ( -  c  ntJS  I  PH )  *S  IN  (  PH2  )  )-  ( COS  ( PH  )*  ( IS  IN  I  P.H3 ) )  **3  U-  CPI -PH  )  /  5125 

1 .  0  )  513.0 

I  5 ( I )  =  XN5 ( I  )  5135 

7  0  n-  I?I1N  1  I  NUF  5140 

UPTURN  5145 

run  5150 


! 


4 

] 


os)  - ; 


(!•/ 


%  Frikf RAN  DECK  7001 

fiPHNil  PUNCHES  THE  A  1C  MATRIX  IN  FORTRAN  FORMAT  SO  THAT  IT  CAN  BE  7002 

C  USED  IN  CUFA  -  FI  UTTER,  ANALYSIS  BY  COLLOCATION  METHOD  7003 

C  THE  AIC  Matrix  FOR  EACH  1/KR  CONTAINS  the  F0LL0WIN6  CARDS  -  7004 

r  f’ARIi  i  -  1/KR  .COLUMNS  1?J.2,  UP1E12.5  3  7006 

n  CARD  2  -  NS17F.  NPART,  NFORM*  NRflW,  COLUMNS  3 -4,5-R>9-12>13-i6  7007 

T  (414),  NFORM=NROW=l  7008 

C  IHF  FOLLOHINO  CARDS  ARE  REPEATED  FOR  EACH  NON-ZERO  PARTITION  IN  7009 

C  IHF  A I C  MATRIX (AS  MANY  TIMES  AS  THE  NUMBER  OF  STRIPS  INTO  WHICH  7010 

r  IHF  SURFACE  IS  » I V I RFll 7011 
r  CAPIi  i  -  Nr  COLUMNS  1-4,  (14)  -  0RDFR  OF  PARTITION  7012 

C  HU  | OHED  HY  1  OR  2  CARDS  FOP  EACH  ROW  OF  PARTITION  MATRIX,  7013 

Cc  rOPMA  T (1 P6E1-?  •  5T  7014 

r  7(iiy 

SllftROHTT NF  PlUiJ  (EKR,Nli,CH,  TSZ,NSEO)  702lj 

!i  IMF  NS  ITlU  Np(T ) »  CH  (4 ,8 , 25  )  7025 

1  FORMAT ( 3  PIE 1 7.5>  6ilX,  3HSTP-,  f  2 , 3HG 0 1  )  7U3U 

?  FOPMAT  (414  .  56X,3HSTP,  I2.3H002)  ‘  703*| 

3  t'nPMAT(I4,fif.y.3HSrP,I2,I2,Il>  7040 

4;  FORMA  r  (iP6r  r?.5i3HSTP,  I?,  I2,Tt)  7044 

5  F0PMATMP2E1?.S,48X,3HSTP,  12,  I?,  M)  7046 

6  F  ilPMA  T-f  1  P3E1  2 .5, 36X,  .TRS^f  P,  12,1?;,  II  )  7047 

7  F OPM  A.  T  (T  P4E1-?  .5, 24X,  3HS  FP» 12, l?i II )  704^ 

0 IJ W H H  .3,'EKR,  NSFQ  7050 

AS  T  2F=0  7055 

•10  H  1=1,  IS/  '  7060 

B  1S.!;?F  =  NS1  ZF  +  N!I(  l )  7065 

•'FOPPrl  7070 

’IRON  =1  7075 

“IIMCH.  2,NSTZF,TSZ>.NFORM,;NROW,NSEO  7080 

•Mi  2  li  1=1.  I 7085 
1SF=3  7090 

>isNII(q)  7095 

:,nnr.H  3,H,NSEIJ,  I.nSF  7100 

‘'SI=NSF  +  1  7105 

7310 

"0  19  K=1,N  4  71,15 

ir(FKim3i14,14  7120 

13  I  F  (N-?)22,21..??  *  73  25 

2,1  r-UNfH  5.  (CH(K,|  ,  J  ),|  =1 ,  H2,?  ) ,  NSEQ,T  i  NSF  7126 

00  10  ih  7129 

2?  I  F ( L-3 )?6, 23 . 26  7130 

??  'MiMCH  6 » ( CH  (>K ,  l  ,  I  ) ,  l  =1 ,  N2 , 2  ) ,  NSFO,  I  * NSE  7131 

.'il  10  IB  73  32 

26  MiMfH  7»  <CH(K,I  ,  I  V,l  =l,N?,7).NStC,  I*NSE  713J 

•  If)  10  18  -  71.34 

14  'FLKEQ.4)  r,n  TO  15  7135 

ir (N-?)25,24 ,?5  713* 

24  ‘MIMCH  7 , .( CH  ( K  ,  I.  »,F  ) ,  I  =]  , M2 ), NSEQ,  I.NSF  713> 

00  TO  18  7331 

28  iMIMCH  4  «  ( CH  <  F ,  I  » I ) , I =1, N2  ) »  NSEQ,  I ,  NSF  714  (i 

00  TO  18  ,  714b 

15  °llMi:H  4 ,  ( C  H  <  K" ,  l  ,  I  ) ,  1  =]  ,.6*)',  NSE  0, I.NSF  71,56 

•iSF  =  NSE  +  l  7355 

PUNCH  5,  (CH(F.I  ,  I  ),l  =/,«),NSE03,  I  ,NSE  73  6(1 

lfi  :SF=NSE+1  716b 

19  ’’ONI  I.NUF  7171 

20  CriMllNUF  7175 

-F  TU^N  71811 

'  END  7l8b 
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I  FORTRAN  13-27-69 

CBJY1  BESSEL  FUNCTIONS  OF  1ST  AND  2ND  KIND,  ORDER  i 

C  FUNCTION  8 Jl,  BY1 

FUNCTION  BJKYY) 
niMENSlON  HES1C5) 

DIMENSION  AC7).BC7>,FC7),FC7>  1 

DATA  ( A ( I >, 1  =  1. 7  )/i 00001109, -.0O03i761».OQ443319,-. 03934289, 

1  .21093573, -.56249986,. 5/ 

•UTA  <B( I ), 1=1,7)/. 0027873, -.1400976,. 3123951, -1.3164827, 

1  2.1682709, .2212091 > -.6366198/ 

DATA  CPU >» 1=1, 7)7. 00029166, -.00079824, -.01074348,. 00637879, 

1  -.00005650!.-. 12499612,. 78539816/ 

DATA  ( F( I ), 1=1, 71/-. 00020 033,. 10113653, -.01249511,. 00017105, 

1  . 01659667,. 156E-5,. 79788456/ 

HAT  A  HES1CD/27H  BY1  ARG.  LESS  THAN  10**-18  / 

X  =  YY 


BJY10I 

BJY1II 

BJYIOO 

BJYlli 

BJYIOO 

BJYiOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYiOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 


IF  (X)  2,1,2 

1  MJ1=0.0 
RETURN 

2  IF  (ABS(X)rJ.n )J, 3,4 
C  BJ1  -  X  LESS  THAN  3.0 

3  7  =  (X#X)/9.0 

30  SUN1  =  All) 

*10  31  1=2,7 

31  SUHla  Z«SUHl  +  ACn 
HJL  =  SUK1*X 
RETURN 

C  BJ1  -  X  GRFATER  THAN  3.0 

4  M»1‘ 

00  TO  42 

C  RY1  -  X  GREATER  THAN  3.0 

41  N=? 

42  ■)  a  ABS(X) 

7  =  3.0/0 
SUM1  =  PCI) 

SUH2*  FC1) 

>10  43  1=2,7 
SUH1=  Z*SUH1»P(I) 

43  SUM?  =  Z*SIIM2»r(!) 

IF  CN.E0.2)  GO  TO  45 

44  SUM  =  SUM2*S  I N  (  U-SUH1 )  /SQKT'C 0  ) 
HJ1  =  SUM 

IF  (  X  •  LT.  0  •  O’)  RJ1=-SUM 
RETURN 

45  RJ1  =-SUH2»CGSCX-SUMl  )/S0RTCX5 
RETURN 

ENTRY  RYl(YY) 

X  =  YY 

6  IF  (X-3.0)  61,61.41 
61  IF  (X-l.E-18)  7,7,611 

C  BY1  -  X  LESS  THAN  3.0, 

6il  /  =  C  X*X,)/9. 0 
SUM1  =  A,C1) 

SUH?=B<1  ) 

9.0  62  1=2,7 
SIJM1  =  Z*SUM1 +AC I  ) 

6?  SUM2  =  Z*SUH2>BCI) 

SUM  =  X*SUM1 

HJ1  =  SUM2/X  +  2.  *Al. OG ( . 5*X  ) «SUM 
RETURN 

7  CAI  l  F  XFM ( 57 , MFS1 , 5 , 5 ) 


BJYIOO 

BJYIOO 

BJYiOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYiOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

HJYlrifl 

BJYIOO 

BJYIOO 

BJYIOO 

BJY100 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

BJYIOO 

1.14159265  BJYIOO 

BJYIOO 
(40>  BJYIOO 


RETURN 

END 


BJY1B861 

BJY10062 


(41) 


f 


*  FORTRAN  83-27-##  RJY008I 

CBJYI1  BESSEL  FUNCTIONS  OF  1ST  ANO  2ND  KINO,  ORDER  •  BJYOOOl 

C  FUNCTION  8 JO,  B  YR  8JTOOOO 

FUNCTION  BJKYYl  BJYOOIO 


niHENSION  HES1(5)  BJYOOIO 

01  MENS  ION  A(7);P<?),Fl7),B(7>  BJYiloi 

PATA  (AM),  1  =  1, 7)/.  0582100, -.8039444, .0444479, -*3163866, 1.2656298,  BJYOOOl 
1  -2.2499997,1./  BJYOOOl 

OAJ A  (PM ), 1=1.7 )/-.lJ558E-3,.29333E-3,.54i25E-3*-. 08262573,  BJYtOoi 

1  .3954E-4, . 04166397; .78539816/  BJYOOOl 

OATA  <F( I), 1-1.7)/. 0001447 6,-.ie0728l5, . 90137237, -.9512E-4,  BJYOOOl 

1  -i 0 055274, -.77E-6, .79788456/  BUYUOOi 


OATA  (B( I ) » 1 =1.7)/-. 00 024846,. 00427916, -.U4261214,. 2530 0117. 
T  -.74350384,,.  60559366,. 36746691/ 

OATA  MFSKl  )/?7H  8tO  ARC.  LESS  THAN  1U*«.-18  / 

X=YY 

IF  (X)  1,3,1 
3-  hjo  =  r.o 

return 

1  IF  (A8S(X)-3.fl)  2,2,4 

C  /X/IS  LESS  THAN  OR  EQUAL  TO  3.0 

2  7=  ( X*X )/9. 0 
5UM=A {  1 ) 

00  21  !=2,7 
21.  SUM=Z*SUH*AH  ) 

RJO=SUM 

return 

C  X  IS  GREATER  THAN  3.0 
4  -)=ABS(X) 

7=3. 0/(5 
SUM1*  P(l) 

SUM2*  FIX) 

•TO  41  1=2,7 
SIJM1  *(  Z»SUM1*P(  \ ) 

4  t  SUM?*  7»SUM2*r<I) 

SUM=SUM2«(:oS(U-SUhl)/SQHT(Q) 

•IJflsSUM 

RFTURN 

a  by  n 

ENTRY  BYO(YX) 

<  =  YY 

IF  (X-3.0)  7,7,6 

6  7  =  3..0/X 
SUMlsp(l) 

SUM?*F(r) 

00  61  1=2,7 
SUM1=  Z*SUM1  +  P(  |  ) 

61  SUM?  =  Z*SUM?  +  MD 

RJO  =  SUM2*STN( X-SUM1 )/SORT(X> 

RETURN 

7  IF  (X-1.E-18)  5,5,8 

8  SUM1 =B ( 1 ) 

SUM2=A(1 ) 

/=  (X#X)/9.0 
00  81  1=2,7 
SUM1  =  Z*SUM1*B( I ) 

8,1  SUM2  =  Z*SUM?*A(!) 

RJO  =  SUN1  ♦  ?.;*ALOO(  .5*X)#SUM2/3.1415926S 
RETURN 

5  CAM  F X E M  ( 5 7  ,  ME S’l ,  5 , 5  ) 
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3.0  PISTON  THEORY  AERODYNAMICS  PROGRAM 
3.1  THEORETICAL  DEVELOPMENT 

The  'pressure  on  s  lifting,  surface  is  normally  given  by  a  surface  func- 

2- 

tional  relationship.  However,-  in  the  limits  of  high  Mach  number  (M  >>  1):  or 
high  reduced  frequency  (M^k  »  1  or  M^2  >>  1) ,  this  relationship  becomes  a 
point  function.  As  a.  bohsequence  of  this  limit,  aerodynamic  influence  coef¬ 
ficients  (AICs)  may  be  Specified  exactly  by  a  strip  theory,,  so  only  a  single 
strip  need  be  considered  in  the  basic  development,  and  control  surface  and 
camber  effects  may  be  determined  in  2  straightf orward  mahher. 

The  present  formulation  derives  the  AICs  from  third-order  piston  theory, 
for  a  (paraholically)  cambering  aiTfoil  with  or  without  a  (rigid  chord)  con¬ 
trol  surface.  The  derivation  differs  only  slightly  from  that  of  Ashley  and 
1 

Zartarian  in  that  in  the  present  case  the  third-order  pressure  coefricient 

is  generalised  to  account,  for  sweep  and  steady  angle  of  attack,  and,  follow- 

2 

ing  a  suggestion  of  Morgan,  Huckel,  and  Runyan,  a  correction  (optional)  is 

suggested  to  give  agreement  with  the  second-order  quasi-steady  supersonic 
3 

theory  of  Van  Dyke.  This  quasi-steady  correction  should  extend  the  validity 
of  piston  theory  to  lower  supersonic  Mach  numbers,  at  low  reduced  frequencies. 
The  derivation  given  here  is  a  combination  and  generalization  of  those  given 
in  Ref,  4  for  the  rigid  chord  airfoil  with  control  surface  and  in  Ref.  5  for 
the  paraholically  cambering  airfoil  without  control  surface; 

The  . AICs  are  defined  to  relate  the  surface  to  the  aerodynamic  forces  at 
the  same  points  in  the  oscillatory  case  by 

{F}.  «  p»2b2s[Ch]{h} 


and  in  the  steady  case  by 


M  ■  &»«>  [°h.]  IM 


Development  of  the  General  Oscillatory.  Case  Including  Camber  and  Control 


Surface 


We  wish  to  determine  the  AICa  that  relate  the  four  def lections  h^,  h^. 
h^,  and  to  the  forcee  F^,  F^,  F^»  end  F^  acting  at  the  same  point*  aa  shown 
in  Fig..  1.  If  the  airfoil  has  no  control  surface  x^  may  be  located  arbitrarily; 
if  there,  is  a  control  surface  x^  oust  be  located  at  the  hinge  line  and  h^  is 
the  deflection  of  the  trailing  edge.  The  deflections  are  those  of  the  mean 


camber  line  and  the  control  surface  is  assumed  to  be  rigid. 


We  consider  two  airfoil  cross-sections.  The  first  is  typical  of  air¬ 
foils  employed  in  missile  applications  while  the  second  is  representative  of 
aircraft  applications.  The  first  airfoil  consists  of  three  straight  lines  as 
shown  in  Fig.  2.  The  second  consists  of  two  tangent  parabolas  and  a  straight 
line  as  shown  in  Fig.  3. 
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y-tg.  ^ .  1 .2  —  First  Airfoil  Idealization 


Fig.  3.1.3  -  Second  Airfoil  Idealization 
The  equations  for  the  deflection  curve  of  either  airfoil  are 
h  »  [h(x)/h1]h1+  [h(x>/hj,]h2  +  [h(x)/h3]h3 

for  0  £  x  £  x^,  and  on  the  control  surface 

hc  -  [hc(x)/h3]h3  +  [hc(x)/h4]h4 

for  x^  ^  ic  Sc  where 

<h(x)/hx  -  (xrx2)(x-x3)/(x1-x2),(x1-x3) 

h(x)  /h2  -  (x-x^  (x-x3)  / (x2-x1).(x2-x3) 

h(x) /h3  -  (x-x1) (x-x2) / (x3-x1) (x3-x2) 

hc  (x)  /h3  »  (c-x)/  (c-x^) 


■j 

i 


4 

3 
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and 


b.c(x)/h4  «  (x-^/Cc-Xb) 

The  equation*  for  the  semi-thickness  distribution  of  the  first  airfoil  are 
,g1(x)/c  «  (T/2)(x/xm)  ,  0SxSxfl 

g2(x)/c  -  (t/2)[1  -  (l-rh)(x-xm)/(xh-xm)]  »  \  \ 

and  g3(x)/c  *  (rh/2)[l  -  .(l-rt)(x^>/.(cij^)]  ,  J^Sx^c 

where  the  symbols  ere  defined  in  the  Nomenclature.  The,  semi-thickness  equa¬ 
tions  for  the  second  airfoil  are 

gL(x)/c  *  (T/2)(x/xm)(2-x/xm)  ,  O  ixix^ 

g2(x)/c  -  (t/2) {i  -  (1  -rh) [ (x-xm) / (xh-xffl)  j2}  ,  x^  *  x  * 

and  g3(x)/c  »  (rh/2)[l  r  (l-rt)(x-xh)/(c-xh)]  ,  ^  Sc 

The  linearized  lifting  pressure  coefficient  for  small  disturbances  about 
the  trim  angle  of  attack  «0  is  given  in  Ref.  4  by 

Cp  -  -  (4v/M)[C^  +  20^  +  3C3M2(gJ  +  a2)] 
where  the  coefficients  ,  C2>  and  C3  are  discussed  in  Ref.  4  and, are /given  by 
\  -  M/(M2  -  sec2A)^/2 

C2  -  [M^(yhl)  -  4- sec2A(M2Tsec2A)]/4\i<2-sec2A)2 
C3  =  (yhD/12. 


and  the  dimensionless-  harmonic  disturbance  downwaah  v  is 


If  the  secant  of  the  sweep  angle  A  of  the  leading  edge  is  taken  as  zero  the 
usualpiston  theory  is  obtained;  if  sec  A  is  taken  as  unity  no  sweep  cor¬ 
rection's  made  to  the  quasi-steady  supersonic  result. 

The  necessary  derivatives  for  the  calculation  of  the  pressure  coef¬ 
ficient  are  the  following. 


Sx  *-h  +  Ch *00 /H23h2 

+  [h “’(x) /h^lhj  ,  0*  x  £  Xj 

and 

H  *  Chc'(x)  /h3]h3(  +  [hc'(x)./h4]h4  x3  ^  x  i  c 

where 

h'(x)/h1  -  (2x-x2-x3)  /(Xj^-x^  (xx-x3) 

h  (x)/h2  *  (2x-x1rx3)/(x2-x1)(x2-:i3) 
h  (x)./h3  «<  (2xrx1-x2)/(x3-x1)(x3rx2)' 

hc'(x)/h3  »  -l/Cc-x^) 

and  hj(x)/h4  -  l/Cc-x^ 

For.  the  first  airfoil  the  thickness  derivatives  are 
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1 


*lx/c  ■  <T/S«>  • 


Q  i  x  £  x 
m 

®2x/c  "  -<T/2)a-rh)/<VV  >  V  *  *  *  *h 

s3x/c  *  ’(Th^2)C1‘rP^c"xh^  »  \  *  x  s  c 

and  for  the  second  airfoil  they  are 

g1v./c  *  (t/x  )(irx/x  )  ,  O^x  ^x 

ix  m  m  m 

8^/c-  -Td-^Hx-xJ/C^)2  ,  Xi*xSxh 

g3x/c  "  "(Th/2)  (:1-“rt)/(c-Jch)  »  Xb  *  x  *  c 


Then  the  pressure  coefficients  in  the  three  airfoil,  regions  are 


Cpl ■  -  -(4v/M)[C1+  2cpglx  +  30^(8^  +  «2)'] 

Pp2  “  "(4v/M)  [Ci  +  2C2Mg?x  +  3C3M2(g^  +  a2) ] 

Cp3  “  -(4vc/M)[C.  +  2C2Mg3x<+  3CyM2(g2x  +  «2)]; 

From  the  principle  of  virtual  work  as  applied  in  Ref.  5  and  from  the  defini¬ 
tion-  of,  the  AICs  the  control  point  forces  are 

Fi  -  qAy{Jom  Cpl[h(x)  /h^  dx  +  Jxh  C  2[h(x)/hi]  dx 

m 

c 

+  Cp3[hc(x) /h^  dx}  ,  i  -  l/,2,3,4< 

5=1 

Therefore,  the  elements  pf  the  fourth  order  AIC  matrix  elements  for  the 
strip  are 

(Ch)ij  =  -(2/M)(l/kr)(Ay/s)|(C1  +  3C3H2e2)(jXh  [h(x)/hi][v(x)/hj]  dx 

+  J  xhthc  (x)  /hi^vc /h  j  ^  dx) 

/  x  x 

+  2C2m|J  “  gjx[h(x) /h-i][v(x)/hj]  dx  +  J  h  g2x[h(x) /h^][v(x) /h^]  dx 

C  Xm 

+  S  83}^[h  (x)/h,][v  (x)/h  ]  dx\ 
xh  '  3  ' 


X  X‘ 

J  ”  gjx[h(x)  /hi'][v(x)7h .]  dx  +  j  h  g2x[h(xi/h  J[v(x)/h  ]  dx 
,o  J  x  J 


+  J  g3x  [hc  (x)  /hjj  [v,c  (x)  /h j  ]  dxj  , 
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where  we  note  that  i,j  «  1,2,3,  and  4  and  also  that  h(x)/h^  «  0  for  i  ■  4,. 


h-(x)/hi  -  0  for  i  «  1  and  2.  We  define  the  following  definite  Integra: 


J1 


-  i(k/b)I^-(§,TD  -  J  g^(x)[h(x)/hi][v(x)/hj]  dx 


so  that 


R^(§,T))  «  j  g^(x)[h(x)/hi][h '(x)/hj]  ’dx 


Jl 


and 


“  /_  8^ (x)  [h (x)  Au] [h (x)  /hj]  dx 


Then  the  AIC  becomes 


^h'lj  '  •<2/MHl/^)(Ay/*)|(C14-3C3M2^)iR[°)(0,xm) 

+  l(k/b)I^CO,y  +  R^vV  +  «R/b)lJ°?(xo,xh)^ 

+  Ecij(\-C)  +  t(k/b)I*°j(i!h,c)]  +  2C2M[R<2)(0,xm)  +  l(k/b)I™  (0  .xj 

*  *U  <W  +  ‘WlgV/k'  :|; 

i 

+Rcu<xh-c>  +  : 

+  3C3M2[R<^(0,*b)  +  l(K/b)l'2)(0,xm)  +  *®<vV  +  i(k/»I<f  Oc^)  • 


+  R$<V'>  +  t(k/b)l'fJ(xh.,=)]| 


where  i*j  «  1,2,3,  and  4,  and  we  note-  that  »  I  »  0  if'  I  or  j  *  4 ;  and 

R  ,  =  I  .  ,  *=  0  if  i  or  j  *  1  or  2. 
cij  cij  J 


\  ' 

1; 


1 '  j.i  ^  x - /*i.; 


-  ‘-V  >;a. 
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3.2  PROGRAM  DESCRIPTION' 

A  general  program  to  calculate  a  set  of  aerodynamic  influence  coefficients 
using  piston  theory  has  been  developed;  The  method  is  applicable  to  wings  of 
moderate  to  high  aspect  ratio  and  speeds  in  the  supersonic  regime.  The  analysis 
can  be  performed  for  wings  with  a  rigid  chord  or  a  flexible  chord.  The  effects 
of  a  flexible  chord  are  accounted  for  by  the  introduction  of  parabolic  cambering. 
Parabolic  camber  is  induced  if  a.  bending  mode  is  parabolic  and  if  a  torsion 
mode  is  linear  in  the  region  surrounding  the  strip  under  consideration.  The 
analysis*  cam.be,  performed  with  or  without  a  control  surface.  The  steady  state 
case  is  available  as  a  limiting  case  of  the  oscillating  case  for  use  in  static 
aeroeiastic  analysis.  The  AICs  relate  the  aerodynam ’ c  forces  to  the  surface 
•deflections  through*  the  following  definitions.  In  the  oscillatory  case, 

iFl  •  /"  Pv\2*  [Ch]  {h} 

and  in  .the  steady  case, 

:  *  ®  <s/S)  [Che]  {hj 

The  AICs  are  derived  for  each  strip  considering  the  airfoil  to  have  up 
to  four  degrees  of  freedom:  pitching,  plunging,  cambering,  and  control  surface 
rotation.  The  program  -provides  the  AICs.  in  printed  and  optional  punched-card 
output  format.  The,  punched-card  output  format  is  identical  to  that  required 
as  input  into  the  COFA,  Collocation  Flutter  Analysis  Program  (Ref.l).  The 
program  capacity  is  25.  surface  strips  and  15  variations  of  Mach  Number.  There 
cam  be  as  many  as  20  reduced  velocities  for  each  Mach  Number. 

3.2.1  PROCESSING.  INFORMATION 

A;  OPERATION 

Standard  FORTRAN  IV  processor  system. 

Operable  on  the  GE  635  computer. 

B. v  COp  STORAGE. 

The  program  STRIP  requires  a- minimum  of  20,000  memovy  units  for 
executuib; 

C.  ADDITIONAL  MACHINE  COMPONENTS 

Standard  FORTRAN  input  tape\>(5) 

Standard  FORTRAN  output  print  tape  (6}, 

Standard.  FORTRAN,, output  punch  tape. 
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3.3  PROGRAM  INPUT  INSTRUCTIONS 

UNITS  Since  ail  of  the  input  dimensions  are  geometrical  and  the;  aerodynamic 
Matrix  is  dimensionless,  only  a  consistent  set  of  length  units  is  necessary  -  inches 
or  feet.. 

DATA  DECK  SETUP 

1.  Title  Card  1 

2. '  Title  Card  2 

3.  NVAN,  NCAM,  NFOIL,  NALPHA,  NTAUS 

4.  ISZ,  MSZ,  NOPUNJ,  JSIZE 

5.  sec  x»  br,  s,,  S,  c 

6.  Ay,  b,  f3>  £oi  e*ch  *triP 

7.  T»Tu»-Tt  for  each  strip 

8.  Mach  Number  Senes 
9.,  Alpha  Series 

10.  1/kr  Series 

Item  1  Title  Card  (Any  alphanumeric  character  Col»anns  1-80) 

Item  2  Title  Card  (Any  alphanumeric  character  Columns  1-80) 

Item  3  Control  Card  (Format  1814) 

Field 
Name 
Column 

NVAN  ■  0  No  Van  Dyke  Correctipn  Included 

*  1  Van  Dyke  Correction  Factor 

NCAM  *  0  No  Camber  Effects  Included 

*  1  Camber  Effects  Included 

NFOIL  «  1  Airfoil  1  (See  Figure  3.1.2) 
v  2  Airfoil  2  (See  Figure  3 . 1 . 3) 

NALPHA  *  1  Angle  of  Attack,.  Alpha  is  constant  for  each  strip 
*•  ISZ  Angle  of  Attack,  Alpha  varies  with  each  strip 

NTAUS  *  1.  Airfoil  Thickness  is  identical  for  each  strip 
■  ISZ  Airfoil  Thickness  varies  for  each  strip 
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‘•su-  . 


Field 

Name 

Columir 


Field 

Name 

Column 


Item  4  Control"  Card  (Format  1814) 


1 

2 

3 

4 

5 

ISZ 

MSZ 

N0P.UNJ 

JSIZE^ 

JSIZE2 

1-4 

5-3 

9-12 

13-16 

ISZ  *  Number  of  strips ;<_ '25 

'MSZ  *  Number  of  Mach  Numbers 15 

NOPUNJ  «  0  Output  Punched 

1-  No  Output  Punched 

JSIZE^  ■  Number  of  Reduced  Velocities,  (V/bu)  ,  <2,0 

for  each  Mach  Number  "i"  where  i  *  1  to  MSZ; 


Item  5  Data  Card  (Format  6E12.8) 


3 

'  1 

< 

2 

3 

,  4 

) 

5 

SECLAM 

BR  ' 

S 

CAPS 

CBAR 

1-12 

' 13-24 

*25-36 

,37-48 

49-60 

SECLAM  =  Secant!  where  Xis  the  leading  edge  sweepback  angle 
BR  =  Reference  Semi-Chord 

S  =  Reference  Semi-Span 

-  -  ISZ 

CAPS  =  Total  Planform. Area,  CAPS  =  2Ayb 

1, 

CBAR  =  Mean  Aerodynamic  Chord 
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DELTAY  «  Ay,  Strip  Width 
b  *  B  ~  Strip  Semi-Chord 

"  Z1 ~  Fraction  of  Chord  to  First  Control  Point 

^2  ■  Z2  ~  Fraction  of  Chord  to  Second  Control  Point  (Use  negative  if  there 
ia  a  control  surface  on  this  particular  strip) 

f,  M  Z3  Fraction  of  Chord  to  Third  Control  Point  (Use  hinge  line  coordi¬ 
nate  if  there  is  g.  control  surface  on  this  particular  strip) 

f  '■  ZM  ~  Fraction  of  chord  at  maximum  thickness  point  on  the=  airfoil 

ttl  ** 
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Item  7  Data  Card  (Format  6E12.8) 


Field 

1 

2 

3 

4 

5 

6 

Name 

taul 

TAUHX 

TAUT1 

TAU  2 

TAUH  2 

TAUT- 

' 

Column 

1-12 

13-24 

< 

25-36 

37-48  ; 

4?-60: 

( 

:  61-72 

TAU  “  Maximum  Airfoil  Thickness  Divided  By  the  Local  Strip 
Chord  Length 

TAUH  ~  Airfoil  Thickness  at  Control  Surface  Hinge  Divided  by 
the  Local  Strip  Chord  Length;  if  there  is  no  control 
surface  set  TAUH  *0.0 

TAUT  ~  Airfoil  Thickness  at  Trailing  Edge  Divided  by  the 
Local  Strip  Chord  Length. 


Repeat  for  each  strip  consecutively;-  two  strips  per  card;  continue  on 
successive  .cards  as  necessary. 


Item  8  Data. Card  (Format  6E12.8) 


Field 

Name 


Column 


1 

,  2 

3 

4 

.5 

6 

EMACH x  • 

EMACH2 

t  #  * 

•  «•  * 

^Hmsz 

>  ( 

1-12  ‘ 

;  13-24 

25-36 

37-48 

49-60 

61-72 

EMACH  ~  Mach  Number 


Continue  on  next  card  as  necessary 
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Item  9  Data  Card  (Format  6E12.3) 

Repeat  this  card  or  Series  of  Cards  for  each  Mach  Nucfcer  when 
NALPHA  *  1 


Column  I  1-12  1  13-24  125-36  37-48  49-60  61-72 


ALPHA  *  Angle  of  attack,  a,  ALPHA  is  constant  for  each  strip. 
Repeat  this  card  MSZ  times. 


When.  NALPHA  *  ISZ 


column  1-12 


13r24  1  25-36  | 37-48  I  49-60 


ALPHA  *  Angle  of  attack,  a.  ALPHA  varies  for  each  strip. 

Continue  ALPHA^  on  next  card  if  necessary.  Repeat 
this  card  or  series  of  cards  MSZ  times. 

(Start  new  card  for  each  Mach  Humber) 


Item  IQ  Data  Card  CForaat  6E12 . 8} 


Field 


Name 


Column 


ERR.  =  Reduced  Velocity  Series;  continue  on  next  card  if 
1  necessary.  Repeat  the  above  card  o t  series  of  cards 
MSZ  times. 

(Start  new  card  for  each  Mach  Number) 


3.4  SAMPLE  PROBLEM 

As  an  example  problem,  the  supersonic  AIC’s  are  calculated  at  Mach 
No.  =*2.5  for  the  high  aspect  ratio  swept  back  wing,  shown  below.  The  wing  is 
analyzed,  for  airfoil  nunber  one  with  parabolic  cambering.  The  analvsis  is 
performed  for  the  reduced  frequencies  (1/k)  of  0.0,  2.0,  and  5,0s  A  1/k  «  0 
calculates  the  aerodynamics  associated  with  stady  state  flight.  A  control 
surface  exists  on  strips  3  and  4. 


PROGRAM  INPUT  DATA 


cosA  *  0 

br  =  4.5  ft 

s  =  12.0ft 

S  -  96  ft2 

c  ■  5.0  ft 

14c  »  5.0,  2.0,  0.0 

M  -  2.5 

?1>  4  2’  ^3  be  any  arbitrary  position  when  no  control  surface  is  present;, 

when  a  control  surface  is  present,  and  x,  ^  may  be  arbitrarily  located  and 
must  be  located  at  the  hinge  line.  Iri  both  cases,  however,  and 

should  be  distributed  across,  the  chord-  so  that  the  clamber  can  be  properly 
defined,  e.g.,  ^  =  .20,  t2~  -S0*  and  ? 3  =  -80., 

NOTE:  Negative  "trips  1  and:  4  indicates  a  control  surface  on  these  strips 
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3.5  PPOGRAM  LISTING  _ 

S  FORTRAN  DECK 

CMfllN  PROGRAM  PISTON  -  AERGDYNAHIC  INFLUENCE  COEFFICIENTS  BY 

C  PISTON  THEORY 

C  - !-TH  OR  WITHOUT  CAMBER 

C  J I TH  OR  WITHOUT  A  CONTROL  SURFACE 

C  '!CAH  =  0>  CAMBER  NOT  CONSIDERED  NCAM  =  1,  CAMBER  INCLUDED 
C  'iFOIL  =  1  OR  2  INDICATES  WHICH  AIRFOIL  IS  TO  BE  CONSIDERED.. 

C  'VAN.  IS  THE  CONTROL  FOR  THE  VAN  DYKE  CORRECTION  OPTION  FOR  USE 

c  iith  Lower  supersonic,  mach  nos.  at  low  reduced  frequencies. 

C  ‘iVAN  =  0,  OPTION  NOT  EXECUTED  NV’AN  =  1,  OPTION  EXECUTED 

C  71, Z2  AND  Z3  ARE  PERCENT  CHORDS  OF  1ST. 2ND  AND  3RD  CONTROL  POINTS 
C-  IN  EACH  STRIP. 

C  '  NEGATIVE  Z2  INDICATES  PRESENCE  OF  A  CONTROL  SURFACE  ON  THE  STRIP;,. 

C  THE  3RD  C.P,  MUST  BE  AT  HINGE  LINE  IF  THERE  IS  k  CONTROL  SURFACE. 

C 

"I  MENS  I  ON  TITLE(24)»DELTAY(25  );,  B  (25),>TAU(25),Z1(25),  Z2 ( 25 ) »  Z3  (.25 ) , 

1  72A ( 25  >  i,TAUI:(  25  ) ,  TADT(  25 ).»  X  ( 25) i  xi(25),X2(25),X3(25),XM(25  )., 
2'M(25;)C,  JSIZF  (15  )  ,EMACH ( 15 )  » NU ( 25  ) 

I  MENS  I  ON  EKR( 28.15),  ALPHA (25.15 ) , CHI ( 4 , 8 , 25 ) , 8H2 ( 3, 6, 25 ) , CH3 ( 3, 6, 
1’5  )  i- CH4  (2 . 4.25).  TN^(  3, 2  ) .  RHO  ( A .  4  ) ,  RATH .  4  ) ,  RA2  < <4. 4  > .  RBI  (4,4), 

2 H 0  ( 4 , 4  )  ?/SAl  (4,4),  SA2  (  4, 4.) ,  S81  ( 4 , 4  ) ,  SB2  { 4 »  4  ) ,  RCO  ( 4,4) ,  RCl  (4,4); 

£  .JC2  (4,4  )/» SCO  (4,4),  SCI  (.4,4-) ,  SC2  (4,4),  T(>4, 3 ) ,  RB2 (4,4  )  ,  D1 ( 3 , 8  ) , 

4  <2(2,6)>TT('3,4),T N'T( 2,3) 

r 

'MJMMON  X, X1,X2,X3,XM„TAU, TAUH.TAU I, RHO, KA1.RA2, kB1,HB2,SH0,SA1, 
1-SA2,  SB1 ,  SB2 »  RCQ ,  RCl ,  RC.2 ,  SCO ,  SCI ,  SC2 
r. 

1  FORMA T( 181 4T 
2,  »0RMAT(6E12,8) 

3  f ORMAT (.iHO  26X, 66HAER0DYNAM TC  INFLUENCE  COEFFICIENTS  BY  PISTON  THE 
1  <RY  WITHOUT  CAMBER  ) 

4  lUPMAfdH  1-7  X,  H5H  (THE  VAN  DYKE  QUASI-STEADY  THEORY  IS-  USED  TO  DET 
l'RHINE  JHE  AERO'iY'JAMlC  COEFFICIENTS.)) 

5  FORMAT  (IHO  54X,TQHINPUT  DATA  //1H  44X,  12,  7H  S1R.IPS/1H  44X.I2, 

11 3H  HACH  NUfi.TFRS  /I47.28H  REDUCED  VELOCITIES  (40IAL)  /1HU  45X, 

2 1 5HSE.CANT  LAMHD*  ,=  E14.6  739X,  22HKEFERENCE  SEMI-CHORD  =  £14.6  /III 
3  49X.11HSEMI-SPAN  -  Ei4.6  /i  H'  46X,=1 4HSURF^ACE  AREA  =  E14.6  /1H  53 

4 v, 7HC  BAH  =  <R14.6  /-IHO  4X.9KSIRIP  NO.  ,10X,7HDELTA  Y , I?X, 1HB, 17X , 

5  MIZl,L6Xi2HZ2,T'>X,2HZ3,i4X,.4H2HAX//(  19,  E24 . 6, 5E18 .6  ) ) 

A  i  ORMATT IHO  42X,9HSTRT?  NO. ,10 X, 2-0 M ALPHA  ZERO  (DEGREES )//( 46X,  12, 

1.’  4 X,  FI 0 . 2  )  ) 

8  FORMAT (1HU  48X, L3HMACH  NUMBER  =  F14.6//(1H  53X,8H1/K(R)  =  E14.6)) 

23  FORMAT (12A6). 

24  MIRMAKIHB  29  X ,  63H.AER0D  YN  AM  I C  INFLUENCE  COEFFICIENTS  BY  PISTON  THE 
1  -RY  WITH  CAhBFR) 

25  ‘•'ORmATUHO  5  J  X',  J  6'H  OSCILLATORY  CASE/-/, XU  45X,li)HMACH  NU .  =  F14.6, 

1//1H-  47X,8»t’/K(R)  =  E14.fi  //I57.7H  ST  RTFS1) 

27  CORHATCIHI  T2A6//1X, 1 2A6/// ) 

28  f ORMATT1HO  4*,9HSTR1P  NO. , 12X, 3HTAU, 13X, 6HT AUT H ) , 12X, 6HTAU ( T )// 

1( 19, E 24.6, 2E18.fi )) 

29  FORMAT (1HU.  49.X,3HCH(  12, HH)  SI7E  =  I2>,3H  BY  12  //) 

30  F  OPMAT  (1H  ?E16 . 8 ,  IX ,  2E16 . 8 , 1 X ,  2E.16 . 8 ,  IX ,  2E 16 . 8  )• 

31  r  ORMAT  C 1 Hi > 

32  FORMAT ( 1P1E)  2 .5 , 6il X,  3HPTN,  I2.3HII01) 

33  FORMAT  ( 4  1 4 ,  *'6X,  3HUTN,  I2.3H002) 

34  FORMA T ( 1 4 , 68 X , 3HP TN , 12, 3  2, II) 

35  FORMAT MP6E12. 5, 3.MPTN,  12,  12,11) 

36  forma  TT 1  P2E.12 . 5 , 48X ,  3HP T.N ,  I2,SI2,  I  I  > 

37  FORMAT (1P3EJ  2.5,3AX^3HPTN,  I2,4'2»  1.1) 

38  FORMAT ( 1 P4E12 . 5, 24 X, 3 HPT N, 1 2 »  12,  II ) 
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11 

12 

13 

14 
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18 

19 

20 

21 

22 
23 
29 
30. 

o 

35 

40 

5U 

55 

60- 

61 

65 

66 

70 

71 

72 
75 
74' 

75 

76 

77 

78 
80 

85 

86 
8  7 
80 
9U 

91 

92 

93 

94 

95 

96 
9/ 

98 

99 
10  U 
lul 
102 


(63) 


39  FORMAT { 1H0  53X,  UNSTEADY  CASE  //1H  45X»10HMACH  NO.  =  F14.6, 

1//1H  47X,  17H1-/K  (  R  )  =  INFINITY  //I57~,7H  STRIPS) 

4  0  C()RMA  K  1 H  29Xi4El0.8) 

41  pORHATUHO  42X»44HT HICKNESS  INTEGRALS  CALCULATED  FOR  AIRFOIL  ii 

42  FORMA  I ( 1HO  42X, 44HTH I CKNESS  INTEGRALS'  CALCULATED-  FOR  A JRFQIL  2) 

•FAD  INPUT  I«ATA  AND  PRINT 
200  ^E AI)  ( 5  >  23  )  (  T  I  TLE < I) , I =1 .24 ) 

•>r  Af)(5,  l )  NVAN,NCAM,NFOILiNALPHA/NTAUS 

'•h  AIK  5,1)  ISZ,M$Z,NOPUNj, ( JSIZE(K),P=1,MSZ) 

>6*11(5,2  )  SFCI.A*i» UR,  S,  CAPS.CBAR 

EAIK5.2)  ( DELTA  Y(  i ) ,  H  (  I ) ,  Z1  C'l  ) ,  Z2  ( T  ) ,  Z3,(  I )  >  ZH (I)  ,T  =1;  I  SZ> 

>fcAD { 5,2)  (  TAUU  >,1AUH( I )iJAUT(  U,  i=l»NTA«JS) 

•  KIT  E  ( 6, 27  )  ( T  LTLELi ) ,  I<-1*  24 ) 

‘ADDEG  =  3.14159265/1.80.. 

i  F ( NTAUS-1 )  224,224,226 
?24  •(]  225  I  =  1  ,<l  SZ 

'  AIK  I  )  =  TAU(1 ) 

1  AtlfK  D-TAUU  LI ) 

?25  1  A.HT (  I >)  =  TAU1  ( 1°).< 

226  'FAIX5,2)(ELACHU;).,  I=1,HSZ) 

•II  227  H=1,.“SZ 

227  >EAIK5,2)  ( ALPHA ( 1 , M ) ; I =1 , NALPHA ) 

( F  (  NCAH.EC).  )r)  GO  10  250 

RITE  (6,3  )< 

:()  TO  251 
250  --RITF  (6,24) 

«5T  f  r(MVAN;EO.I!)  GO  TO  252 
■RI1E(6,4) 

V52  i  F  (  NF 0  I.L.  EO  •  2  )  GO  TO  253 
iR  f  IF  (6,41=)- 
1(1  254 

?53  R IT  F ( 6 , 42 ) 

75  A  iSIIM  =  0 

■  0  255  I=l,i-S7 
755  i SUM  =  JSIJM  ♦  JSIZE(I) 

"0  300  1=1, 1 SZ 
^ 2 A ( I )  =  ABS ( Z2 ( 1  )  ) 

> (  I  )=2 . 0*B (  I  ) 

■'i  ( n=x  ( i  )»zi  ( i  > 

2(1  )  =  X(‘I  )*Z2A(  I ) 

>i( F)=X( I )»Z3< I ) 

« 0  0  vH(  |  )  =  X(  I  )#/.M(  I  ) 

RI  1  E  (  6, 5  )  ISZ.-lSZ,  JSUM.SECL  AM,  UR,  S,  CAPS.CBAR,  (  I ,  DELTA  Y  (  I  ) ,  B  (  1)  , 
171(1  ),  Z2A(M,  Z3(  I  ),ZM(  I  ),  1=1,  ISZ) 

RITE (6, 28  )  ( I , TnU( I )  ,  rAUK( I ), TAUT (  I ), 1=1, ISZ) 
iF(NALPHA-1  )236, 236, 2.58 
:36  -0  237  1=1, ISZ 
■0  237  M=T,  f.SZ 
237  »LPHA( I,H)sALPHA(1,M) 

738  U  240  1=1, "SZ 
iS7  =  JSIZE( I ) 

■>F  AO  (  5 , 2  )  (EKR<  J,  I  ),  J=l,  JSZ) 

R  IT  F  (  6 , 8  )  IMaC'I  (  I  ),  (EKRU,  I  ),  J=1,JSZ) 

R I J  E( 6, 6 ) ( J, AL^HA ( J> I ), J=l, ISZ) 

•  0  240  J=l, ISZ 

240  ALPHA!  J,  I  )=  ALPHA  (  J,  1  )*RAI)UEG 
0  1000  M=1 , MS/ 

.*MS  =  EMACH(H)*I  h/\CH(M;) 

^f:cs=sfclam«sfci.am 


116 

117 

118 
120 
121 

129 

130 

131 

132 

133 

134 
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140 
14.1 

142 

143 

144 
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147 

148 

149 

150 

151 

152 
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155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 
167 
16H 
170 
171 
25  0 
255 
260 


C  '‘AN  DYKE  OPTION 

IF(NVAN)  310,310,312 
<10  «:HAR1=1. 0 

-:BAK2  =  (3.4*1-fd)/4.0 
,  :()  TO  320 

.<12  'HAHl  =  EHACH(H)/SOKT.(EHS«SECS.) 

:BAR2=7<EHS*I  HS*<2.4)-4.fl*SECS*(EHS-SECS> >/ ( 4 . Q* ( EMS-SECS ) » f EHS- 
l'ECS)) 

320  *:RAK3=2.4/1?. 

• SFU=0 

•S7= JS ! ZE ( H ) 

•n  900  J-l/JSZ 
rF(tKR(J,M) 3325,330,325 
'25  (1  -=  1  .0/EKK(J,.'<) 
r?si..n/(Fl*H  ) 

■' 3  =  Fl/BR  T 
130  0  FOO  1=1,  IS7. 

C  ^TFAOY  CASE  OPTION 

iF  ( I- KR(\J,M)  3340,335,340 
^35  ■f'ON  =  (4./EHACH(H)  )* (CBAR*DLLT AY  (.IT/CAPS) 

-  Ml  10  344 

<40  UN  (<2./EKACH(M)  )*F2*<DELTAY<  IT/S) 
r.  I  RF  OIL  OPTION 

344.  IF  { NFO I L  «  EG»2 )  GO'  TO  350 
IF  (Z2(  U.BT.O.t.  »  00  TO  345 
'AIL  TH I N1  { 1,J  ) 

:t)  TO  360 

<4**  AIL  THINK  0,1) 

-0  TO  360 

<50  i  F  {  Z2  ( 1  )  •  iLT»  0 . 0  )  00  TO  355 
ALL  THI N? ( 1  ,IJ) 

0  fO  360 

355  '.ALL  TKP’\2  .(Oil  ) 

<60  (I  50  0  K  =  l,4 

0  400  .L  =  l,  7,2 

IF  (  K  .  E G  ♦.!  #  UR  «.K  •  E  0  •  2  )  00  TO  361 

■HI.  TO  362 

■361  •  F(l.  .EQ.7  3  DO  To  399 
•50  TO  365 

<6?  -F  ( K . E  0 • 4  )  GO  TO  363 
:0  TO  365 

163  IF  (  L .  EO  •  1 »  OR .  1.  •  E  D  •  3  )  00  JO  39y  * 

;()  TO  365 

<99  :H1  vK,l,  I  )  =  i'.U 
50  TO  400 
<65  i.L=(L  +  l)/2 

0  -•'ASIC  A I C  HATIM-a  tUOAl'ION  FOR  P1SI0N  THEORY  (CAMBER  WITH  A  CONTROL 
C  UURTACF)  -  REAL  ELEMENTS, 

CHI  (  K ,  I. ,  I  )  =  -CUN*(  ( CB AR1  +  3 # *CBAR3* 

1PMACH(H)*«2«AI  PhA(  I,M)*«2)#(RH0(K,LL)-»RC0(K,LL)  )  +  2,*CBAR2*EMACH(M) 
2*  (RA1(K,LLMRH1  (K,LL.)  +  RC1  (K,LL  ))  +  3,*CBAR3«EMACH(H)**2*(RA2(K*LL)  + 
3-’H2  ( K,  LL )'4 RC2  ( K , LL  ) )  ) 

<1(10  :ONT  I  NIIF 

D  450  L=2 ,1-,? 
i  F  (  EK R  (  J ,  M )  )S7 0 , \i\ 9 , 3  70 
370  IF  (  K .  E  0 « 1 » (‘  R «  K .  E  M  »  2  )  GO  TO  371 
.0  TO  372 

<71  >.F(L  ,E0*8)GO  TO  409 
•JO.  TO  375 

372  ( F ( K , E0.4  )  (-0  T'l  373 


264 

265 
270 
275 
2B0 
285 
2  9  0 
291 

305 

306 
3io 
315 

319 

320 
325 

33  0 
335 

34  0 
345 
350 
355 
360 

404 

405 
4  06 

409 

410 

411 
415 
4.19 
42.U 

421 

422 
425 

435 

436 
4  37 
43H 
43y 

440 

441 

442 
44  3 

444 

445 

446 
446 

449 

450 
455 
460 
465 
4  711 
475 
4  HO 
481 
4  82 
4  83 
484 
4  85 

(65) 


V it  - 


*0  7(1  375  4Q6 

.<73  'F<L.E0w2.0K.L.6Q.4>  GO  TO  404  •  487 

*•0  7  0  375  488 

4»9  ‘HI (  K , L , I  1  =  0 . 0  489 

:0  10  450  490 

<75  • L=L/2  491 

r  :ASIC  AiC  MATRIX  EOUAIIGN  FOR  PISTON  THEORY  (CAHBER  WITH  A  CONTROL  493 

C  ‘aiRFACE)  -  IMAGINARY  ELEMENTS.  494 

CHl(K.L.I)  =  “CON* ( (CB AH1+3 • *CBAR3#  495 

1FMACH(H1**2*ALPHA(I,H1»*2)«F3*(SHU(K>LL)+SC0{K»LL))+2.*CBAR2*  50  0 

2f  MACH  (  H  ).*F3*  (  SA1  ( K »  LL  )  +  SBl <  K >  LL  )  +  SC1  { K *  LL,) )  +3 . *CBAR3*EHACH( K ) **2  505 

3iF3«(SA?(K,lL>+SB2<K,LL)*SC2(K,LL>)>  510 

a  5  0  ‘ONIINIIE  ’  515 

‘•0  0  \OK\JI  NllE  520 

C  GENERATE  AH.  MATRICES  525 

C  CHI,  -  CaMHER  H 1 7H  A  CONTROL  SURFACE  530 

r  CH2  -  RIGID  CHORD  WITH  A  CONTROL  SURFACE-  535 

r  CH3  -  CAHuEB  WITHOUT  A  CONTROL  >SUR'F4CE~  540 

C  CII8  -  RIGID  CHORD  WITHOUT  A  CONTROL  SURFACE  545 

if  (?2(  I  HGT.O.blGU  TO  566  635 

l  F ( N'CA’H.E 6 .-!•')  G!1  10  550  636 

II (  l  }  =4  637 

.0  10  800  639 

‘•51]'  0  660  KK  =  1  j  A  64,0 

0  560  JJ=1,3  645 

iv6  0  ’  (KK,  JJ)  =  0.0  650 

*  C i > 1  )  =1 •  655 

■(?,li  =  (Z2AU)-/.3(I.))/(Z3(I)-23(I))  660 

( ? ,  2 • )  =  (  Z 1  (  !  ) -Z2 A  ( I>) )  /  ( 2 1  ( I )  - Z 3 1 1 ) )  665 

•(3,2)=1.  670 

’ ( 4  »  3 ) =1 «  675 

ALL  MULT  ( 1 . CH  l ,  CH2 »  U1 »  T  T , 3 , 6 , 4»  8 ,  I  )  695 

11(11=3  696 

(1  TO  80  0  698 

•i66  1 0  567  K  =  l-i  t  649 

»0  567  L  =  l,6  700 

r,67  >:h.1(K»L«  I)a|;Hl(4*L»  I  >  701 

T  F  ('NCAM.EO.  ‘I-JfiO  TO  600  702 

■HU)  =3  70  3, 

■0  TO  800  704 

*nn  TN(1'»1)=1.  .  705 

’  N (  1,2  )  =0.0  706 

1  N ( 2»  t  )  =  ( Z2 A  ( I  )-Z3 ( I  /c)  / ( Z1  ( I  )-23 (1)1  •  710 

•N(?,2)  =  (Z1M  )  -  Z2  A  (  I  ’)  )/(Zl(  I  )--Z3  ( I  )  )  715 

rN'(3, 1  )  =  0 .0  7,5  0 

<  N ( 3, 2 1 =1 .  -  725 

ALL  MULT  (  I  N »  C-13 »  C|i4  ,  P? ,  TNT  ,  2 , 4 , 3 , 6  ,L)  745; 

U ( I ) =2  746 

‘’UO  "ONTINIJE  750 

t:  ’R  I  NT  A IC  MATRICES  800 

■  R II  F  ( 6f»3 1 )  805 

i  F  (  NC AMv>  E'J .  1  >  G  l  TO  825  111  II 

H  I  I  E  (  6  »  3  )  li  t  I 

0  TO  830  812 

<25  R  I  7  F  '(  6 , 2  4  )  81.) 

•'30  IF  (NVAN.Eu.ii  )  RT  10  835  814 

RITE(6*4)  ■  815 

•’35  IF (FKR(J.M) 1836,837,836  810 

H36  K I TE ( 6 » 25 1  FMACH ( M ) , EKR ( J, M ) , I SZ  817 

•I)  10  838  818 


(66) 


*37  *RITE<6,39)  EHACH<H>vISZ  819 

H38  HO  675  !=r,JSZ  825 

«  =  NU<I)  826 

'?  =  2*N  827 

•RITE (6, 29)  I,N,N  828 

•0  874  K=1,N  829 

iF(EKR<J,M) )840»8/l,840  .  830 

-MO  iF(Z2(  I  ).GT.O.O)  GO  TO  866  83i 

I F ( NCATt.EQ.  U  ).  GO  TO  850  832 

•  R!.TE(6»30)  (CH1TK,L,  I  ),L=i,.N2)  835 

•0  TO  874  838 

*50  Rl  TE(6,3Q  )  (CH2(K,Li  I  >,L=1,N2)  843= 

';0  10  874  .  845 

H66  JFFNCAM.EQ.U )  G'j  TO  870  848 

RITE(6,30)  {CH3(K»L,  I ) , L=l, N2  )  85,3 

0  10  874  855 

870  .R11EC6,30)  (CH4<K>L» 1 )>.»L=1,N2T  860 

•0  TO  874  ‘  861 

<’71  iFC72(  I  ).GT.jl.,0  )  00  TO  873  862 

irOiCAH.EO.U-i  G«I  ro  872  863 

•HI  TE(6,40  )  (CHl  (  K  >  L »  l  ) »  L=1  , N2, 2  >  864 

0  10  H74  865 

•'■72  RI1EO6.4U0  ( CH2 ( K»  L ,  I ) ,  L=1  ,  N2,  2)  866 

•  *.()  TO  874  867 

*73  i F{ NCAH.EO. 0 )  G«  10  855  •  368 

•Rl  I£(6,4(L)  (CH3(K,L,  I  ),L=1,N2,2)  869 

iO  T 0:  874  H70 

>•55  -R II E ( 6 * 411 )  (CH4(K,L>I>,L=1,N2,2>  871 

<v74  -OMIINUF  872 

*75  ’ONllNUt  873 

C  MINCH  A1G  MATRICES  IN  FORTRAN  FORMAT  SO  THAI  IT  CAN  BE  USED  IN  875 

n  COFA  -  FLUTTER  ANALYSIS  BY  COLLOCATION  METHOD.  878 

C  fHF  A I C  MATRIX  FOR  EACH  1/KR  CONTAINS  THE  FOLLOWING  CARDS  r  880 

r.  f'ARI)  1  -  1/FR,  COLUMNS  i-12»  (1P1E12..5)  883 

C  CARO  2  -  NS  I  ZE>  NPART,.  NFORM,  NRQR,  (414),  NFCRH=NRUk=1  885 

C  IMF  F0LL0RIM5  CaROS  ARE  REPEATED  FOR  EACH  NONrZEKO  PARTITION  IN=  888 

C  TilF  A I C  MATRIX  (AS.  MANY  T4M£S  AS  THE  NUK8ER  OF  STRIPS  INTO  WHICH  89U 

0  THE  SURFACE  IS  ill  V ))f»ED ) .  /  '  893 

C  CAR 0  1  •?  N,  ORDER  OF  PARTITION*  COLUMNS  1-4,  C  14)  895 

C  COL I  OWED  BY  1  OR  2  CARDS  FOR  EACH  ROW  <0F  PARTITION  MATRIX  89/ 

C  MINCHED  FORMAT  <lP6tl2.5>  •  898 

G  899 

•  F  (  NOPUN J  ) 9(i  11 , 8 76 , 9li 0  900 

*76  iSFO=NS£Q*l  905 

=MJNCH  32,EKM  J*M),NSEG  91 U 

SI/E  =  0  '  915 

0  877  1=  1,  ISZ  920 

>  77  :SIZE  =  NS I'ZE  +  Nu(  V)  925 

FORM  =  1<  930 

ROW  =  1  935 

"UNCH  33,  NS  I  ZE,  .1 SZ »  NF URPif.NR.GW,  NSE<J  940 

•0  89,5  1=1,  ISZ  ‘  94  5 

SF  =  1  96(1 

a  THUD  955 

oIJNGH  34,M,NSFQ„I  ,  NSE  960 

SE  =  NSE  ♦  1  965 

2  =  N#2  970 

0  694  K  =  1 ,  f'  975 

■  F ( FKR ( J, M ) )  87  S* 881,878  976 

878  LF(  Z2GT  ).GTiO.-0  >  CO  TO  886  980 


(67) 


r 


IF(NCAH.EQ.O)  GO  TO  880 
PUNCH  35,  (CHI  (K,  L»  I 1,  L=1 , 6 ) , NSEG,  I,  NSE 
SE  =  NSE  ♦  1 

PUNCH  36,  (CH1(K,L,  1),L*7,8),NSEQ, I,NSE 
GO  TO  893 

R80  PUNCH  35, (CH2(K,L, I >,L=i,N2)*NSEG, I,NSE 
.  00  TO  893 

886  IFTMCAM.EO.II )  GO  TO  890 

PUNCH  35,  (CH3CK,L„I  > , L=1 , N2 ) , NSEQ,  f  ,'NSE 


I 


l 


¥ 


: 


% 


L 

L 


L 


oil  TO  893 

H90  PUNCH  38,  (CM  (K,L,I. ), L=1  , N2 ) , NSEQ,  I , NSE. 

50  TO  894 

*>81  I F ( 22  (  I  ),GT>. 0.0  )  GO  TO  883 
>f ( NCAM.EO. b  )  GO  TO  882 

•MINCH  38,  (  Clii  (K,L,  I  ).,  L=1  *  N2 , 2 ) ,  NSEQ,  I ,  NSE 
JO  TO  893 

5*82  JUNCH  37,  (CI.2(K,L,  I  ) ,  L-l,  N2 , 2  ) ,  NSEC,  I ,  NSE 
;0  10  893 

!-83  IF(HCAH.EO.I')  GO  TO  884 

;*UNCH  37i  (Ci«3(K,L,  I  ) ,  L=l,  N2 , 2  )  i  NSEQ,  I ,  NSE 
Ml  10  893 

H84  “UNCH  36,  (Ch4(K,L,  I  ),L=1,N2,2).,NSEQ,  I,NSE 
*<93  SE=NSE+1 

894  CONTINUE 

895  TUMINUF 
*>on  -OUT  I NUE 

l.UO  MONTI a  Lie 

:0  TO  200 
CNN 

%  F OUTRAN  DECK 

CT-IN1  COMPUTES  I  HE  TH I  CIS  NESS  INTEGRALS,  BOTH  REAL  AND  IMAGINARY 

C  FOR  AIRFOIL  lx. 

C 

SURROUT I NE  I H I N l ( NCO, I ) 

T I  HENS  I  ON  X125),XL(25),X2(25),X3(25),XH(25),TAU(25),TAUH(25), 

1  >  A1 (4,4),RA2(4,4), RIM  i4,4),SH0(4,4),SAl(4,4),SA2(4,4) ,  S&l (4 , 4 ) , 

2 MR? ( 4 , 4  ) , RA ( 4 , 4 , 3 ) , RB ( 4 , 4 , 3 ) i SA ( 4 , 4 , 3 ) , S8 ( 4 , 4 , 3 ) , RC (4,4,3), 
3CC(4,4,3),RH0(4,4),  T  AUT<(  25  ) ,  RB2  ( 4, 4  ) ,  RCO  ( 4 , 4  ) ,  RC1  ( 4, 4  ) ,  RC2  ( 4, 4  )  , 
4  :C 0  (  4 , 4  ) ,  SC  1  ( 4 , 4  }/,  SC 2  (  4 , 4  ) 

♦OMMON  X,X1.,X2,  <3,XM,TAU,  I AUH , TAUT , RHO , RAX, RA2, RBI, RB2, SHO ,SA1, 
1MA2,SH1,SB2.RC0,RC1,RC2,SC0,SC1,SC2 


i 


L 

I 

i 

1 

I 


5 

6 


10 


15 


•:l  =  (  TAU(  Ic)/2.  )f  (X(  I  )/XM(  I  ).) 
iFLNCO.EQ.1)  GO  fu  5 
<H  =  X (  l  ) 

i AUH(  L)  =  TAUT(  I  ) 

:0  TO  6 
vH  =  X3(  I  ) 

;3  =-(TAUH(  I  )/2.  )  *  ( 1  •  -TAUT  (  D/TAUH(  I  ) ) *X  ( I  j /  ( X (  IT “XH ) 
;2  =-(TAU  (!)/?. )*(1.-TAUH(I)/TAU(I))«X(I)/(XH-XM(I)) 
■>U  100  11=1,4 
iO  ion  j=i , 4 

Mj  T0(10, 15, 20,50),  I  I 
I  =  X  i  (  I  ) 


••'K  =  X2 ( I  ) 
■>'l  =  X3  (  I  ) 
'! 0  TO  30 
v I  =  X2 (  I  ) 
/K  =  X1(  I  ) 
v  L  =  X  3  (  I  ) 


985 

990 

991 

992 
995 

1000 

1005 

1010 

1015 

1020 

1025 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1038 

1039 
104  0 

rb4i 

1042 

1U43 

1044 

1045 
1050 
1055 
20  OU 
2001 
2002 
2003 
2005 
2010 
2011 
2012 

2013 

2014 

2015 
2016. 
202,0 
2025 

2030 

2031 
2035 
20  4  0* 

2045 

2046 
2.055 
2065 
2070 
2075 
2U6U 
2  0  65 
2090 
20  95 
210  0 
2105 
2110 


(.68) 


,0  TO  30  2115 

20  1 1  =  X3 ( I )  2120 

-K= XI M  )  212b 

'  L=X2  (T).  ,  2130 

30  0  70(35,40,45,50), J  2190 

35  <J=Xl(f)  2195 

VP=X2(I)  2200 

*0=X3(I)  '  2205 

■:0  TO  55  221U 

4  0  vj=X2(I)  2215 

*P=XiH)  2220 

•0=X3(i>  2225 

0  TO  55  .  2230 

4*5  <J=X3  ( I  )  2235 

<P=X1.(J)  2240 

■'0=X2U)  2245 

-  0  10  55  2250 

50  'HfM.I  I,.J)=0.0f  '  2255 

•A1  ( I  l,  J)=0.0  2260 

:'A2  (I  L, J )  =  0 . 0  2265 

HI  (  II,  Jl  =  £i.n  2270 

B2 ( I i , J ) =0 .0  2275 

tAl  (  i  I,  J)=6.0  22BO 

:A2  C I  L,  J ).f  0 . 0  2285 

*.RJ  ( IT,  J)  =  0.0  •  ,2290 

"H2 ( I ! , J ) =0 . 0  .  2295 

•HO  (!  T,  J)  =  0.0-  2300 

'F(HCO.fcCI.O)  *30  TO  90  2301 

;F( II. £0.3. OR. II.EO/4)  UO  TO  54  2302 

:0  10  90  2305 

54  1  F ( J. EO. 3. OP.J; EG. 4 )  GO  TO  100  2306 

'  :0  TO  90-  '  2307 

55  El  =  <  XI  —  XK  )  *  6X I*^XL  )* <  XJ-XP )* ( X J«*XG  )  2310 

0  75  K=l,3  2315 

0  TO  (60, 62, 64), K  2320 

60  ■!  =  1.0  2321 

J2  =  1  .0  2322 

:0  TO  69  ,  2324 

62  *1  =  03  2325 

•>2  =  G2  2326 

•  0  10  69  2327 

64  -1  i  61*61  2329 

■'?  =  62*62  2330 

69  :«F(K-1) 71, 70,71  ’  ,233o 

70  X  =XH  233/ 

TO  72  2330 

71  'MAX  =  X M (  I  )  23,39 

72  A(  II,J,K)  =  (P1/,JEL><*  (  X M A X* #4/ 2  «  -1  *  /3  .,« ( 2  c  *  XK +2  »'*XL+ XP ♦  X Q  )  * XM A X  2340 

1**3*. 5  *i2.*XK*XLKXK+XL)*(XP  +  X0)  )*>MAX  **2-XK*XL*(  XP+XO  )*XHAX  )  2345 

*.a( n,d, k)=(p.!/med*( (Umax  #*5)/5. )-(xk*xl+xp+xo)*( (xhax  **4)/4.)  2350 

1  «  FXK*XL  +  XP*xb+( XK  +  XL )*( XP  +  XO ) )*< ( XHAX  **3)/3. )- ( XK»XL* ( XP+XU ) +  2355 

2  P*X0*(XK  +  XL)  )*(  (  X  M  A  X  *  *  2 )  /  2..T  ♦  X  K  *  X  l  #  X  P  *  X  Q  *  XM  A  X  )  2360 

‘  F  (  K- 1°)  73 , 75 , 73  2361 

73  ■»(  I  J',K),=  (P2/IEL)*(  (XH**4-XM(  I  )  **4)  72  . -1 . /3  .  *  ( 2.  *XK  +  2 .  *XL  +  XP*X(1 )  2365 

1  *(Xh4*3-XM'.(  I  )*#'.>’)  *  ( 2o»*XK*xL+  (  X  L  +  X  K)  *  t>X  P  +  X  Q  )  )  *  (  XH**2-XM  ( l>)*«2)/2.-  2370 

2  'K»XL*(XP  +  XM).*(XH-XM(  l  )h  ”  2375 

.114  I  I,.JVK)  =  (P2AHEI.)*(  (XH»#5-XM(  I  ).«#5  )/5i-(  Xk*XL*XP*XQ)*<  XH#*4-  2380 

1 -HI  I ) **4 ) /4 >  +  ( XK* XL*  XP*X0* ( XK+XL ) * ( XP+XU ) ) * ( XH**3-XM( i>*»3)/3.-  2385 

2<  XK*Xl.*(XPf  YQ)  +  xp*Xi,*  (  XK* XL  ).;)*(  XH«*2-XM(  I )**2 )/2i ♦XK#Xl*XP*XQ*  2390 

3 ' XM-XM ( I ) ) )  2395 


(69)- 


1 


75  -nNlIMUO  2400 

-A1 ( II, J)=RA( 1 1, J,2)  2405 

•A?<  PI,  J)=RA(  IJ*J,3>  2410 

’fM  ( 1 1 » J)=R8(  1 1*  J,2)  2415 

R?(II»J)sRHCM,J>,£)  2420 

A1  (  I  li  J)  =  SA(  1 1,  j,2)  2425 

a?( IF, J)  =  SA( 1  I, J,3)  2430 

•'HI  {  1 1,  J)=SH<  II„  J,2)  2435 

H2  (  I  I  ,  J)=.SH  (  II , J, 3 )  2440 

■'HO  ( 1  I ,  J )  =  RA:(  1  I  , J, 1 )  2445 

:HO(I  I:,  J)=SA(!I»  J,l)  2450 

•T'{  l-I  .E0.3.ANI).  J.EO,.3)GO  TO  79  2455 

it  TO  90  2459 

79  .f(NCO.EQ.l)  ISO  TO  80  246U 

:(J  10  90  2461 

Hfi  (l  65  KK=ti3  2465 

•  n  T  0 1 81 ,  83  •  BP.)  >  KK  2470 

81  *3  =  1.(1  ’  2475 

:il  10  84  2480 

83  ’3  =  153  2481 

0  TO  84  2462 

HP  -3  =  63*63  2485 

84  -■C(3,3#KK)^-P>3/2.  .2.490 

*C  ( 4  ,>3 » K  K  )  =  I*  C  ( 3 » 3 »  K  K  )  2495 

C  ( 4 , 4  ,  KK>)  =-RC  (  t,3»KKJ  •  2500 

T ( 3 , 4 »  KK  )  =9(J (4 , 4 , KK )  2505 

C(3,3,KK)=l;3*(  MI  )-XH)/3.  2510 

.C ( 4 , 4  , KK ) =SC ( 3, 3 *  KK )  2515 

6(3,4, KK)=hl*< U I )-XH)/6i  2520 

C ( 4 ,3 , K  K ) =SC (3 ,  4,‘Kk )  .  2525 

85  HPIINUF  2526 

•0  68  LH=3,4  2330 

0  H8  LN=3,4  2535 

CO(LM,LN)=l>C(L'«U  N,l)  2540 

•i:j  ( £’  M >  L N  )  =66  (  L  *l,  i.N*  2  )■  2=545 

CPU  M,LN)=PCU.o,lN,3)  255U 

•Cfi(LM,LN)=SC(L'i,LN,l}  2555 

:C1  (LM,I  N)  =  MC(L*5.»LN,2)  P560 

rC?<LM,LN)=f,C<L'1,LN,3)  2565 

68  '’ON'f  I  NjJt  2566 

■(I  TO  10  0  257  0 

9!)  ‘COM  I,  J)  =  0.0  2575 

;C1  (  11,  J)=0.0  2.68U 

CP  C  I  1.8  J)  =  0.0  2685 

C0(II,J>=0.0  2690 

C1(M,J)=n.n  2-595 

C? (II , j)  =  0 • 0  2  b  0  u 

00  'ON  TIME  •  2605 

fc  T  URN  261 U 

•  np  2615 

*  FORTRAN  PECK  3UUU 

nT','iN2  COMPUTES  I  HE  THICKNESS  INTEGRALS,  &0Tri  REAL  AND  IMAGINARY/  3001 

r  fltR  AiRfOIL  2.  3  H  H  P 

0  30  03 

•HlRKOllT  I  Nt  IHIN2(nC0,  J»)  3  ll  06 

TMF  NSION  X  (  25  )  i  XI  ( ?5  1/  X2  (25  )  /  X3  (  25'j  />XM  ( 25  )  /  TAU  ( 25  )  /  TAUN'l  25  )  /  3Ulu 

1  HO  (4, 4)/ RAJ  (4„  I),Ra2(4,4  ) ,  RB1<  4 , 4  ) ,  Rb2  <  4 , 4  > ,  SHU  <  4 >  4  ) ,  S A!1  (  4 , 4  ) ,  3  01 J 

2  A?  (  4 , 4  ) ,  SB1  (  4 ,  '») »  SH2  (  4  #>4  ) ,  KCO  { 4 , 4  )  /  RC1  ( 4 , 4  ) ,  KC2  (4  /  4  ) ,  SCO  (.4/4),  3  111  P 

3  .Cl  (  4 , 4 ) ,  SCP  ( 'i ,  4  )  i  RC:<;4  >  4  /  3  )  /  SC  ( 4  ,  4/  3  ) ,  X  J1  ( 4 , 4 , 2  )  ,  XK1  (  4  ,  2  )  /  3  013 

4  J?  (  4 , 4  /  2i)  ,  x KP  ( .  \ ,  4 , 2  ) , R4  <  4 , 4 , 2 )  ,'RH  (4,4, 2  ) ,  S A  (  4 , 4 , 2  > ,  SB  (  4 , 4 , 2  ) ,  3014 


(70) 


5aJ0(4,4),AJ1<4,4  ) , AJ2 (4, 4 ) , B JO <4 , 4 ) , B J1 (4,4 ) , BJ2 (4, 4 ) , AKO ( 4,4) >  3015 

6AK1(4,4)»  AK?  (4,4),BKO(4,4),BKi(4»4},BK2(4/4)»TAUTf25)  301t> 

^OHMON^  X*  XI, X2,  X3»  XM*TAU,. TAUH,  TAUT,  RHO,  Ra£vRA2>  RBI,  RB2, SHO / SA1,  301? 

1<:A2,S81,S82,RC3»RC1,RC2,SC0,SC1,SC2  '  3018 

3019 

IF(NCO.EO.I)  GO  TO  5  3020 

VH=X ( I j  3025 

IAIHK I  r=TAUT( I )  3030 

;0  TO  6  3035 

5  'H=X3(I)  .  ,  3040 

•:3‘  -r  (  TAUH {  1  )  /2.  )*M. -TAUT  ( 1  l/TAUH < IT ) *X (J  >/ ( X( I ) -XH >  3041 

6  '1  =  1AU(1)*(X(I)/XM(I)):  3045 

‘2  =  TAIJ{  i )  *  (1  .  -  TAUH  (  I  )/TAU(  1 3 )  *X  <  I)>XH(  I  )/(XH-XH(  I  )  )**2  3050 

•;H1=-Hl/XH(  I  )  3055 

■H2=-U2/XM(J)  3060 

I>  150  1 1=1,,  4  3065 

0  150  J=l,4  30?o 

«l  10(10,15,20,50  ),  I  1  ’  307  5 

in  \!  =  X1:(I)  3080 

*'K  =  X2  { I  )  3085. 

’L=X3(I)  3090 

0  TO  30  •  3095 

15  I =X2 (I)  3100 

•K=Xt(I)  3105 

»L=X3(i)  .  •  3110 

0  10  3o  3115 

20  I =X3  C I )  3120 

•K=X1CI)  3125 

L=X2 (I )  .  3130 

30  :(j  10(35,40,45,50 ),J  .  !  3135 

35  J=X1M)  31 4  0 

P=X2 ( I  )  3145 

*«J=X3(n.  3156 

;0  10  55  3155 

40  *  J=X2 ( I )  3160 

•  P  =  X1 (4  )  3165 

(J=X3  ( I  )  3170 

0  TO  55  3175 

45  »’  J  =  X 3  (  I  )  3180 

'P-Xl(I)  3185 

■  ll=X2  (  I  )  3190 

HI  TO  55  3195 

50  H(l  (  I  I  ,  J  >  =  0  .  fl  3205 

U1(II,J>  =  0.0  321  ft 

-‘A2  ( 1 1 ,  J )  =  1) .  0  321  5 

’HI  CT  I,  J)  =  0.0  3220 

K?(  I  I,  J,)  =  0.0  3225 

■HQ  (  M,  J)  =  0.0  3230 

<A1  (  I  I ,  J  )  =  0 . 6  3235 

•:A?(  1  I,  J  )  =  0  .  0  3240 

:B1(1I,J)=0,0  3245 

■:B?(  I  1,  J)  =  0.(1  3250 

■  MMCO.EO.O)  GO  TO  135  3251 

1  F  (  1  I .  F 0 » 3 .  •> R .  I  1 .  t: 0 . 4  )  GO  TO  54  3252 

.0  10  135  3255 

54  :F(  J.E0.3.0H.  J.t.0.4)  GO  TO  150  325o 

;(1  10  135  3257 

55  F.L=(Xl-XK)*<XI-XL)*(,XJ-XP)*(XJ-XO>  326u 

(1  75  K=l»2  3265 

PU.&Q.2)  i.O  Id  /O  327.U 

(71) 


.'HAX  =  XH  327b 

•0  10  71  3280 

7ri  •  HAX=XH( l >  328P 

71  -A{  U,  J,K)  =  (1./:iEL?*CXMAX*«4/2.-1./3.«{2.*XK+2.«XL+XP+XU}*XHAX»*3  3290 

1  ‘1  i/.2.*(2.«XK»Xl.*(XK*XL)*(XP*XQ)  )*Xr.AX**2-XK*XL»CXP+X0J*XHAX>  329b 

‘;A  (  1 1»  J/I^)  =  H  i/»lEL)*(  ( tXMAX««5)/5.  )-{XK+XL  +  XP  +  X(l)*C(  XHAX#*4)/4i)  330  0 

1*  IXK*XL+XP#xQ*UK  +  XL)*<XP*XQ)  )»(  {XHAX#e3)/3.)-(XK*XL»(XP*XQ>+  330b 

2^P*XQ*(XK  +  X!  )  )*(  (X^AX*»2X/.2.)*XK#XL»XP*XQ*XHAX)  3310 

•F(K.EO.i)  Pfl  TO  75  3355 

•»(  1  l,  J,K)  =  Cl./iJEL)«(  (XH**4-XMAX**4)/2,-il./3.  )*<2.*XK+2.*XL*XP+Xu  3320 

1  •*<XH**3-XH.'X**3}  +  (2.*XK*XL*(XL+XK>*(XP  +  XG)  )* (  XH**2-XMAX**2  )/2. -  33?b 

2  K*XL*(XP  +  X«))*(  xH-XKAX) )  '  3330 

H{  1  I,  J,K)  =  (l./-lEL)*(  CXH**5-XMAX*«5)/5»-{XK>XL+XP*X0)*<XH**4-XHAX  333b 

1**4  )/4.+vXK«XL**P*XU*(Xk*XL)*(XP+XQ) >«(XH**3-XHAX**3)/3,-(XX»XL*  334 11 

?-Xp+Xy)  +  XP*XCl*(  XK*Xl  )  )*{XH**2-XHAX**2  )/2.  +XK*XL:*XP*XO*  ( XH-XHAX ) )  334  b 

7b  0N1INUP  33  b  i) 

M0(  1 1,  J)=RM  I  I.  J,l)  0<^., 

HiH  I  I,  J)=SM  I  I,  J,l)  ‘  ^46Ui 

JO ( 1 1  *  J )  =  Ra  (oj  I ,  j,,  2 )  356!, 

JBC  1 1.  J)=R}«{  { I ,  J ,  2  )  33?  i) 

Kft(  |jj  J)=Sft  (  1 1 ,  J,2)  337., 

KtM  i  I,  |  I,  J>2)  3  38u 

0  »>5  L=l,  2  33hb 

•  F  ( l  .  EQ >2 )  i.O  To  80;  j39ll 

Z=,l?n  •  .  339b 

N  =  XH(I)  34HU 

0  10  340b 

‘  7  =  XH  (  I  )  J4iu 

N=rXH  3  4i  b 

01  *.J1(  l  I »  J/t  J  =  (l./nfcL)*(  .4#{'XN«*5-X2**5)-.25«(2.*XK+2,*XL+XP  +  XUJ*(  34'2,U 

1  ■  N *  *  4  -  X  l *  *' *>  +  <  ?„  *VK.*„X L  ♦  <  X  A  ♦  X  L  )  *  (  X  P  +  X  Q  )  )  *'(  X  N  *  *  3  r  X  Z  *  *  3  )  /  3 ,  -  ,5  *  X  K  #  X  L  *  34  2  b 

2«iXP  +  XU)*(XN»*?^xZ**2)  )  '  0430 

K1  (  I  1>J  >;L  )  =  (l./iDhL)*l{XN*«6-XZ**6)/6.-;2*(XK  +  XL'‘-XP  +  XU)*(XN**b-  34  3*> 
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